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Abstract 


The understanding of the physical processes that lead to the origin of matter in 
the early Universe, creating both an excess of matter over anti-matter and a dark 
matter abundance that survived until the present, is one of the most fascinating 
challenges in modern science. The problem cannot be addressed within our current 
description of fundamental physics and, therefore, it currently provides a very strong 
evidence of new physics. Solutions can either reside in a modification of the standard 
model of elementary particle physics or in a modification of the way we describe 
gravity, based on general relativity, or at the interface of both. We will mainly 
discuss the first class of solutions. Traditionally, models that separately explain 
either the matter-antimatter asymmetry of the Universe or dark matter have been 
proposed. However, in the last years there has also been an accreted interest and 
intense activity on scenarios able to provide a unified picture of the origin of matter 
in the early universe. In this review we discuss some of the main ideas emphasising 
primarily those models that have more chances to be experimentally tested during 
next years. Moreover, after a general discussion, we will focus on extensions of the 
standard model that can also address neutrino masses and mixing. Since this is 
currently the only evidence of physics beyond the standard model coming directly 
from particle physics experiments, it is then reasonable that such extensions might 


also provide a solution to the problem of the origin of matter in the universe. 
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1 Introduction 


The discovery of the Higgs boson at the Large Hadron Collider (LHC) |1| has confirmed 
the Brout-Englert-Higgs mechanism |2| for the origin of masses of elementary particles in 
the standard model of particle physics and fundamental interactions (SM), providing the 
last missing piece for a full confirmation of its validity. Since then, signs of new physics 
have been eagerly awaited at the LHC, in particular those that could have finally provided 
(directly or indirectly) indications on the nature of the dark matter in the Universe, but 
so far no compelling evidence has been found]! 

However, despite null results at colliders, the existence of new physics is strongly 
supported by the necessity to reconcile the SM and the cosmological observations within 
a unified picture. In addition to dark matter, another strong indication of new physics is 
given by the failure to find, within the SM, an explanation of the absence of primordial 
anti-matter in our observable universe, clearly indicating CP violation to a macroscopic 
level. 

These two cosmological puzzles have been traditionally addressed separately but clearly 
a unified picture for the origin of matter in the universe is not only greatly attractive but 
also quite reasonable. In the last years many ideas have been proposed. In this review we 
will primarily concentrate on extensions of the SM also able to explain neutrino masses and 
mixing, currently the only evidence of new physics coming directly from particle physics. 
Therefore, if the origin of matter in the universe can be regarded as a phenomenological 


guidance toward a more fundamental theory, vice versa neutrino masses and mixing could 


'The most recent measurement of the double ratio of branching fractions in B meson decays, Rg = 
0.84610 -044, by the LHCb experiment exhibits a 3.1o deviation from the SM expected value (Rx = 1), 
hinting to a violation of lepton universality [8]. This can intriguingly be regarded as the indication 
of the existence of new physics at a scale ~ 100 TeV suggesting that the discovery of new particles 
might be within the reach of colliders in a near future. Moreover, the value of the muon anomalous 
magnetic moment recently measured by the Fermilab Muon g-2 Experiment, combined with the previous 
value measured at the Brookhaven National Laboratory E821 experiment, deviates at 4.1o from the SM 
predicted value jal. This is calculated combining perturbative expansion in the fine-structure constant 
a with non-perturbative approach based on dispersive relations that require experimental information 
on the hadronic cross section of e+ + e7 [5]. However, both anomalies require further investigation. In 
particular, the deviation of Rx from the standard model prediction is still not statistically significant 
enough considering that it represents just one anomaly among a multitude of other measurements not 
showing signs of new physics. Moreover, recent lattice calculations find an hadronic contribution that 
results into a predicted value of the muon anomalous magnetic moment in the standard model that is 
in good agreement with the measured value [6]. It is then premature to claim discovery of new physics 


without first a clear understanding of the correct standard model prediction. 


hold the solutions to the cosmological puzzles. 

This is the plan of the review. In Section 2 we briefly summarise the observational 
evidence for the matter-antimatter asymmetry of the universe and why this should be 
regarded as a strong indication of new physics. We also provide a general overview 
of models of baryogenesis. In Section 3 we summarise the phenomenological picture 
that strongly points to the existence of a dark matter component in the universe of 
non-baryonic nature, unaccountable within the SM or even within a minimal extension 
where neutrinos are massive Dirac fermions|?| In Section 4 we provide a brief review of 
neutrino masses and mixing and how this can be described within extensions of the SM 
based on the simplest type-I seesaw mechanism. In Section 5 we review the main results 
on leptogenesis including some recent ones. In Section 6 we show how the role of dark 
matter can be played by a long-lived right-handed (RH) neutrino with keV mass produced 
from the mixing with left-handed (LH) neutrinos within the type-I seesaw mechanism, the 
Dodelson-Widrow mechanism. We also discuss how the active-sterile neutrino mixing can 
be consistently embedded within the type-I seesaw mechanism with neutrino mass and 
mixing experimental results from neutrino oscillation experiments. The picture, dubbed 
as yMSM model, can also potentially address the matter-antimatter asymmetry puzzle 
with leptogenesis from RH neutrino mixing realising a unified picture based on a minimal 
extension of the SM. The vMSM can be tested experimentally in different ways and 
stringent constraints have been found. Currently, it is not clear whether viable solutions 
still exist within the parameter space. In Section 7 we show how the introduction of 
a 5-dim operator where RH neutrinos can couple to the Higgs boson can induce a RH 
neutrino mixing producing an amount of long-lived ‘dark’ RH neutrinos able to reproduce 
the measured dark matter abundance. In addition the matter-antimatter asymmetry can 
be reproduced by two RH neutrino resonant leptogenesis. Finally, in Section 8 we draw 
some conclusions, discussing the future of research of models on the origin of matter in 


the Universe. 


2 Baryogenesis 


If big bang nucleosynthesis represents the first example where knowledge of nuclear physics 
could be applied to cosmology and helped understanding the existence of a hot early stage 
in the history of the universe, as confirmed by the discovery of the cosmic microwave 


background, with baryogenesis there is an important step further: the entire universe is 


?This does not exclude explanation of the origin of matter in the universe where neutrinos are Dirac 


particles. However, adding just Dirac neutrino masses is not a sufficient ingredient of new physics. 
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used as a unique laboratory to test new theories, a way to complement ground laboratories 
and go beyond their limitations. Here we discuss the basic points and review some recent 
progress, mainly aiming at showing how the matter-antimatter asymmetry of the universe 


is an important motivation and investigative tool to go beyond the standard model. 


2.1 Matter-antimatter asymmetry of the universe 


The cosmic baryon abundance is today accurately and precisely measured by an analysis of 
CMB temperature and polarisation anisotropies. In particular, the ratio of the amplitude 
of the first to the second acoustic peak in temperature anisotropies is particularly sensitive 
to the value of the baryon abundance}? The latest analysis of the Planck collaboration, 
combining Planck results on CMB anisotropies and baryon acoustic oscillation data, finds 


for the baryonic contribution to the energy density parameter 
Q goh? = 0.02242 + 0.00014. (1) 
Using the Planck result for the Hubble constant 
h = Ho/(100 km s™* Mpc~') = 0.6766 + 0.0042, (2) 
one has for the baryonic contribution to the energy density parameter 
sgo = 0.049 + 0.001 . (3) 


Acoustic oscillations are not sensitive to the sign of the baryon number so they actually 
measure the total amount of matter and anti-matter rather than the net one. If we 
consider the possibility that at the time of recombination there is also some amount of 
anti-matter in the form of anti-nucleons, CMB would be then sensitive to the sum of the 


two contributions. One can then easily derive the total baryon-to-photon number ratio 


= Q 7 
nify = BOT NE a CED a aa 1071". (4) 
Nyo Mp Nyo 


Notice that in this expression, in general, one might have a spatial dependence of ngo 
and ngo but still a homogeneous sum, as acoustic peaks require. For example, one can 
have matter and antimatter spatial domains in the universe. However, at the moment, let 
us exclude this possibility and consider the homogeneous case. Searches of antimatter 
in cosmic rays, coming from space, find no evidence of primordial antimatter. The 


amount of antimatter that we observe in cosmic rays, mainly positrons and anti-protons, 


3An insightful discussion on the physics of microwave background anisotropies can be found in [r]. 
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can be explained in terms of astrophysical processes. This implies that at the time of 
recombination one can assume the amount of anti-nucleons to be negligible. The observed 
baryon abundance then necessarily corresponds to the net baryon asymmetry that had to 
exist prior to recombination and that survived the consecutive stages of annihilations of 
particles and antiparticles taking place in the early universe while the temperature was 
dropping down. We can then find from Eq. the value of the net baryon-to-photon 
number ratio at the present time, since Nko ~ NBo = (Ngo — Ngo) /Nyo = NBtec, SO that 
Eq. translates into 

ngo = (6.12 + 0.04) x 107". (5) 


This value together with the absence of primordial antimatter in our observable universe 
is in agreement with the well known result [9] that starting with a vanishing asymmetry 
one would eventually be left with equal values of the relic abundances of matter and 
antimatter far below, about ten orders of magnitude, the measured value in Eq. (5). 

However, this argument still does not exclude the possibility, neglected so far, of a 
universe consisting of a patchwork of matter-antimatter domains on scales À = 100 Mpc, 
i.e., larger than those probed by the lack of primordial antimatter in cosmic rays[] If we 
consider such a scenario, implying the existence of some mechanism that could segregate 
matter from anti-matter with the creation of matter-antimatter domains on very large 
scales, one could think of an overall baryon symmetric universe, with the cancellation of 
the contributions from matter and antimatter domains. However, if such large matter 
and anti-matter domains exist, and have comparable geometry, then the annihilation 
radiation would contribute to the cosmic y-ray diffuse background. Since such excess 
is not observed, this excludes the existence of domains of matter and antimatter on 
scales A = 20Mpc, so even larger than those probed by cosmic rays and as large as 
the entire observable universe, thus ruling out the possibility of a baryon symmetric 
universe. Therefore, today, barring some caveats, one arrives to the conclusion that a 
matter-antimatter asymmetry had to exist prior to recombination [11]. 

Primordial nucleosynthesis is also sensitive to the amount of baryonic matter. Within 
a standard picture, with no matter-antimatter domains, the value of the baryon-to-photon 
ratio inferred at the time of nucleosynthesis, when Deuterium forms, is in nice agreement 
with the measurement from CMB. Therefore, the matter-antimatter asymmetry, and in 
particular the baryon asymmetry, needs to be generated prior the onset of nucleosynthesis. 
Possible deviations from the standard big bang nucleosynthesis scenario have to be such 


not to spoil the success of its predictions. For example, one could consider the existence 


4Strong constraints on scales A < 20 Mpc, i.e., within the size of clusters of galaxies, also come from 


y-ray observations of X-ray emitting clusters not showing evidence of annihilations 10). 


T 


of anti-matter domains at the time of nucleosynthesis. However, these, on the basis of the 
arguments discussed above, need to contain a sub-dominant component of anti-matter 
compared to matter and need to be sufficiently small to avoid the constraints on a baryon 
symmetric universe. In this case the synthesis of primordial anti-nuclei is expected inside 
these small anti-matter domains in a similar way to how standard big bang nucleosynthesis 
proceeds in matter domains. Interestingly, the AMS-02 experiment has recently claimed 
the detection of six events compatible with being anti-*He and two events with anti-*He 
nuclei [12]. These events cannot be explained in terms of spallation processes, predicting 
an anti-*He flux two orders of magnitude below AMS-02 sensitivity and a anti-‘He flux 
approximately five orders of magnitude below. If confirmed, this discovery would point to 
the existence of compact antimatter domains in the form of anti-stars or anti-clouds 
[13]. The anti-nuclei should have been created during big-bang nucleosynthesis and 
the measured isotopic ratio, roughly 4He : 3He ~ 1 : 3, would be obtained for 7g = 
np/n, ~ 1-6 x 10-1’. However, the creation of these compact antimatter domains in 
the early universe would be highly non trivial to explain. It would require a very specific 
mechanism within the very early Universe, likely related to the same mechanism that 
generates the asymmetry. For example a modification of the Affleck-Dine mechanism 
has been proposed [15]. We will not consider this kind of models but it should be clear 
that a confirmation of the existence of primordial antimatter would likely rule out the 
mechanisms of baryogenesis we will discuss, or in any case it would require an extension 
able to describe a mechanism of formation of antimatter domains in the universe. 

The measured value of the matter-antimatter asymmetry, in the form of baryon-to-photon 
ratio in Eq. (5), needs to exist prior to the onset of nucleosynthesis but it cannot be 
understood in terms of some value of the baryon charge pre-existing the inflationary 
stage and that remained constant until the present time [16]. This would indeed conflict 
with successful inflation requirements since necessarily the energy density associated to 
the baryon charge had to be higher than the inflation energy density, approximately 
constant during inflation, before some time in past. Before this time inflation would 
be then inhibited, since it requires the inflation energy density to dominate. One can 
easily verify that in this way inflation could not last more than N = 6-7 e-folds and 
not the required N ~ 60 folds needed to solve the horizon problem and to explain the 
amplitude of CMB temperature anisotropies on super-horizon scales. For this reason 
the matter-antimatter asymmetry has to be explained by a mechanism of dynamical 
generation, from processes that occurred at the end or after inflation, a stage commonly 
referred to as baryogenesis [17]. 


2.2 Models of baryogenesis 


There are three necessary conditions that have to be typically satisfied by a model of 


baryogenesis for the generation of a baryon asymmetry at the macroscopic level: 
e Baryon number non-conservation; 
e C and CP violation; 
e Departure from thermal equilibrium. 


These conditions are all satisfied by the original model proposed by Sakharov [17], though 
not explicitly stated there, and for this reason they are commonly referred to as Sakharov’s 
conditions. Although there are many loopholes and there are even examples of models that 
do not respect any of these conditions [18], Sakharov’s conditions hold in most traditional 
models. 

There is a very long list of proposed models of baryogenesis}}| Some of them can be 


grouped within the same class, this is just a partial list (as emphasised by the dots): 


e From heavy particle decays: 


— GUT baryogenesis 28]; 
— Leptogenesis : 


e From a first order phase transition at the electroweak spontaneous symmetry breaking 


(electroweak baryogenesis) ; 
in the SM [80]; 


— in the two Higgs doublet model (31); 
— in the minimal supersymmetric standard model (MSSM) [82]; 


in the next-to-minimal supersymmetric standard model (NMSSM) [33]; 


in the nearly minimal supersymmetric standard model (nMSSM) [84]; 


e Affleck-Dine baryogenesis [14]; 


e Baryogenesis from black hole evaporation 38); 
>There are different excellent specialistic reviews or monographs on baryogenesis : 


e Spontaneous baryogenesiq)| [89]; 
e Gravitational baryogenesis [40]; 
e Cold baryogenesis [41]; 

e... 


Here we will briefly discuss the first two classes of baryogenesis models, those that are 
more easily probed by different experimental tests and, for this reason, have attracted 


more interest in recent years. 


2.3 Electroweak baryogenesis 


Electroweak baryogenesis]'is a popular model, or better a class of models, of baryogenesis 
that was initially proposed within the SM framework. However, experimental findings rule 
out this possibility and, since no other models of baryogenesis relying on known physics 
have been found, the matter-antimatter asymmetry of the universe is currently regarded 
as evidence of new physics. 

At the classical level the baryonic current j is conserved. However, quantum effects 
can give rise to an anomalous current contribution that breaks conservation and 
in the case of the electroweak theory one has 


2 
; g TE 
ô, ji = Nf 3972 G" Guy (6) 


where in the chiral anomaly ns denotes the number of families, Gy is the electroweak 
gauge field strength and G uw its dual. 

The chiral anomaly can be re-expressed as a total derivative „K“ of the anomalous 
current K” that does not vanish at infinite so that integrating from a time t = 0 to some 
final time tp between two vacuum states corresponding to Gay = 0, the change in the 


baryon number is non-vanishing and one obtains 
AB(ts) = ANes = ne [Nos(ts) — Nes(0)], (7) 


where Ncs = f d?x K° is the Chern-Simons number, an integer number. At zero temperature 


the transition rate between two vacua is tiny: [spn ~ e7167?/9 m 107160, However, at high 


Spontaneous baryogenesis is a typical example of a baryogenesis model that does not respect 
Sakharov’s conditions since CP is conserved and the asymmetry is generated in equilibrium. In this 


case the crucial ingredient is a temporary, dynamical violation of CPT invariance. 
TFor specialistic reviews on electroweak baryogenesis see or for a pedagogical introduction 


see [45]. 
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temperatures, corresponding to energies above the energy barrier height separating two 
vacua and given by the static solution of field configuration, transitions over the barrier 
can happen. The lowest energy barrier height is found as a saddle point solution called 
sphaleron [49]. The energy barrier height is found to be a few TeV’s at zero temperature 
and in the case there would not be observable effects. However, at high temperatures 
one has to include finite temperature effects and the energy barrier gets lower until 
it disappears above a critical temperature [80]. During the electroweak phase transition 
the universe would then pass from an unbroken phase, where sphaleron processes are 
unsuppressed, to a broken phase where they are suppressed. The transition rate at finite 
temperatures between two electroweak vacua with AB ¥ 0 can then be expressed as 


aa] 


Tsph x T* exp -« (8) 


where v = (®) is the vacuum expectation value of the Higgs with v = 0 for T = Te 
(unbroken phase) and v = v(Z.) for T < T, (broken phase). 

Therefore, above the critical temperature the baryon number can be violated satisfying 
the first Sakharov condition. If the electroweak phase transition is first order, then during 
the transition broken phase bubbles nucleate. Both inside (broken phase) and outside the 
bubble (unbroken phase) a macroscopic baryon number cannot be generated. The reason 
is that during a sphaleron transition lepton number is also non conserved together with 
the baryon number, in a way that AB = AL = 3. In this way one has A(B — L) = 0, 
while any (B + L) number is washed-out in thermal equilibrium. However, on the bubble 
wall there is a strong departure from thermal equilibrium and if the phase transition is 
strong first order, a baryon asymmetry can survive inside the bubble as we explain in a 
moment. First, we need to discuss how the baryon asymmetry is produced. 

Departure from thermal equilibrium and baryon number violation are still not sufficient 
conditions to create an asymmetry: one also needs that CP is violated, otherwise, on 
average, for any process creating a baryon number there would be a CP conjugated 
inverse process that destroys it. The asymmetry is actually initially produced in the 
form of a chiral asymmetry on the bubble wall since LH and RH particles interact with 
different rate with the bubble wall if CP is violated. In this way LH particles penetrate 
the bubble more efficiently and this generates a chiral asymmetry on the bubble wall. 
It is then crucial that the chiral asymmetry is reprocessed into a baryon asymmetry by 
sphalerons inside the bubble and this implies that electroweak baryogenesis ends when 


the temperature drops below the out-of-equilibrium temperature of sphalerons 


Pow ws 189 Gey . (9) 


sph — 
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It should be emphasised that the sphaleron rate needs to vanish quickly inside the bubble 
once the baryon asymmetry has been produced, since otherwise after having reprocessed 
the chiral asymmetry into a baryon asymmetry, sphalerons would also destroy the baryon 
asymmetry if they are in equilibrium. From Eq. one can see that this requires the 
condition u(T.)/T. 2 implying that the phase transition has to be a strongly first 
order phase transition f] It is quite non-trivial that successful electroweak baryogenesis 
could, in principle, be realised within the SM. The source of CP violation would be in this 
case provided by the presence of phases in the quark mixing matrix, the CKM matrix. 
Moreover the Higgs potential could in principle experience a strongly first order phase 
transition. Therefore, potentially, the observed baryon asymmetry could be explained in 
the SM, since all three Sakharov conditions could be satisfied. However, one finds that the 
value of v(T;)/T, in the SM is related to the Higgs mass and it turns out that a strongly 
first order phase transition can occur only for my < 72 GeV [53], clearly in disagreement 
with the measured value my ~ 125 GeV. 

For this reason successful electroweak baryogenesis necessarily requires an extension 
of the SM. For example, electroweak baryogenesis can be successfully implemented within 
supersymmetric models{™] In the MSSM one obtains an upper bound on the Higgs 
127 GeV and on the stop mass my < 120 GeV [55]. Current experimental 
constraints from the LHC have basically ruled out this scenario now [56], though some 


mass my S 
loopholes might still be possible [57]. Within a non-minimal supersymmetric model, like 
the NMSSM, an allowed region still survives and also in a two Higgs doublet 
model only special solutions are found within the full parameter space [60]. However, the 
recent improved result from the ACME experiment on the electron dipole moment of the 


electron [61], 


ldo] < 1.1 x 107” ecm (90% C.L.), (10) 


essentially rules out also these models (see for a recent discussion). 

These findings of course pose the legitimate question whether electroweak baryogenesis 
is still alive [62]. Introducing a scalar singlet, it is possible to realise easily a strongly first 
order phase transition and at the same time avoid experimental constraints from colliders. 


However, this seems quite an ad hoc model specifically engineered to show that successful 


8Notice that the criterium is approximate and gauge dependent. A more sophisticate and gauge 


independent approach is discussed in |52]. 
Consider indeed that for a second order phase transition one has v(T.) = 0. 
Recently a different strategy based on effective theory has been proposed (54). It has been shown that 


augmenting the SM by higher-dimensional operators, the variation of the SM couplings with temperatures 


is changed and this can lead to successful electroweak baryogenesis. 
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electroweak baryogenesis is still viable. Maybe, if certainly not dead, it is fair to say that 
electroweak baryogenesis is then currently is in a sleeping beauty mode, awaiting to get 


awake from some supporting experimental finding] 


2.4 Baryogenesis from heavy particle decays 


Baryogenesis from heavy particle decays is the oldest class of models of baryogenesis, 
since the original proposal by Sakharov belongs to it. Here we review the main common 
features and then specialise the general discussion to an important specific historical 
example: GUT baryogenesis. In Section 5 we will see that leptogenesis provides another 
remarkable concrete realisation within this general class of models that is tightly related 


to neutrino masses. 


2.4.1 Out-of-equilibrium decays 


Let us consider a self-conjugate heavy (Mx >> Mpw) particle X whose decays are CP 
asymmetric, in such a way that the decaying rate into particles, I yx, is in general different 
from the decaying rate into anti-particles, Tx, and such that the single decay process 
into particles (anti-particles) violate B — L by a quantity Ag-z (—Ag_z). The total CP 


asymmetry is then conveniently defined as 


ix Tx 
EX = sign(As-1) 5 Ty . (11) 
For example, for GUT baryogenesis one has sign(Ag_;,) = +1, while for leptogenesis 
sign(Ag_,) = —1. At tree level ex vanishes. Therefore, a non-vanishing total CP 
asymmetry relies on the interference between tree level and one-loop diagrams, where 
in the loop at least a second massive particle Y needs to be present. This means that in 
the limit where My — 0 one necessarily has £x — 0. Moreover €x also vanishes in the 
case of exact degeneracy, i.e., for My = My. 

The total decay rate T'R =T x +Tx is the product of the total decay width, re, times 
the averaged dilation factor (1/y), explicitly: 


r? =T? (=). (12) 


y 


Tn this respect it is intriguing that a strong first order phase transition might also produce a 
gravitational wave stochastic background testable with gravitational wave interferometers. In the specific 
case of an electroweak baryogenesis strong first order phase transition, the spectrum would peak in the 
millihertz range that will be tested by an interferometer such as LISA (see and references therein). 
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As discussed in 2.3, sphaleron processes, while inter-converting B and L separately, 
preserve B — L and for this reason the kinetic equations are much simpler if the B — L 
evolution is tracked instead of the separate B or L evolution. Moreover, it is convenient 
to use, as an independent variable, the quantity z = Mx/T and to introduce the decay 
factor Dx = T2/(Hz). Another convenient choice is to track the abundance of X 
particles, Ny, and asymmetry, Ng_z, in a portion of comoving volume normalised in a 
way to contain one X particle on average in ultra-relativistic thermal equilibrium (i.e., 
NY(z « 1) = 1) 

The results for the evolution of the X abundance and the B — L asymmetry can be 
conveniently described in terms just of the total decay parameter defined as 

D 
Kx= eo ; (16) 

implying Dx = Kx z (1/7). This is proportional to the ratio of the age of the universe 
when the X particles become non-relativistic, for z = 1, to the X lifetime, Tx = 1/T2. 

The simplest case is realised when Tx is much longer than the age of the Universe 
at z = 1, corresponding to have Kx < 1. In this way the bulk of decays occur when 
the temperature is much below the X mass and the X production from inverse decays, 
or other possible processes, is Boltzmann suppressed. In this situation decays are the 


only relevant processes and the kinetic equations for the X abundance and the B — L 


!2Tn general the abundance of some generic quantity X in a generic portion of comoving volume R? a’ 


is defined as Nx = nx R} a, where Ro is a normalization factor. Our normalisation corresponds to take 


3 _ 4 mT" gs 
° 3 €(3) gx gso 


Another popular normalisation for the abundances is to calculate them in a portion of comoving volume 


(13) 


per unit entropy, assumed to be constant, so that s(T) Ra’ = 1 and Yx (T) = nx(T)/s(T). If entropy 

is indeed conserved, then the two variables are simply proportional to each other, explicitly 
135 ¢(3) gx 

Yx(T) = ——— = Nx(T). 14 

(T) = a Nxt) (14) 

In the calculation of physical observables such as ng (T) the normalisation factor clearly cancels out so that 

all choices are equaivalent. However, the normalisation of the kind NẸ¥ (z < 1) = 1 has a few advantages: 

it does not assume entropy conservation and at T >> mx one has conveniently Ny} (T >> Mr) = 1, while 

_ 135¢(3) gx 


Y (T > Mx) = A gh’ (15) 
S 


not a very convenient value. Notice moreover that often one starts from Yx variables but then 
quantities Yx/Y% (T > Mx) are plotted since they are more convenient. However, clearly one has 


Nx = Yx/Yx(T > Mx) showing that Yx variable is an unnecessary intermediate auxiliary variable. 
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asymmetry are particularly simple to be written, 


dNx 


aX = —Dx(2)Nx(2), (17) 
adNp_p _ dNx 
dz = =EX dz , (18) 
and solved, 
Npg-r(2z) = Ni_p tex [Nx — Nx(z)] (19) 
Nx(z) = Nye aA (20) 


where Ni,_, denotes a possible pre-existing B — L asymmetry and N$ the initial number 


of X particles. The dilation factor, averaged on the Boltzmann statistics, is given by the 


(TRE ey 


and can be simply approximated by 


(asks = 


This simple approximation makes possible to solve analytically the integral in Eq. (20), 


ratio of Bessel functions, 


yielding the result 


z2 fi n pix 
Nee) ove = [2 eee) (23) 


In particular, the final B — L asymmetry is given by 
Nb- = Np_p tex Nx. (24) 


The baryon-to-photon ratio at recombination can then be obtained dividing by the number 
of photons at recombination N** and taking into account that sphalerons will convert 
only a fraction dsp, œ 1/3 of the B — L asymmetry into a baryon asymmetry. In this way 


one can write: ; 
N 
we 2 SBOE 


neg > rec ` 
3 NE 


(25) 


It is useful to introduce the efficiency factor, defined as the ratio of the asymmetry 


produced from X decays to the CP asymmetry: 


Npg-r(2z) 


EX 


Ni_,=0 
K(z; Ky) = BEL 


, (26) 
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— numerical 


analytical 


z=M/T 


Figure 1: Out-of-equilibrium decays for an initial thermal X abundance. 


depending uniquely on Ky. In the case of out-of-equilibrium decays one has simply 
K(z) = NY — Nx(z) and in this way Eq. can be recast as 


Np_1(z) = Ni_,t+ex K(z). (27) 


The final value of the efficiency factor is then simply given by the initial number of X 
particles, i.e., ke = K(oo) = NX. In particular, it is equal to unity in the case of an initial 
thermal X-abundance with z; < 1. In Fig. 1 we show two examples of out of equilibrium 
decays, for K = 107? and K = 1074, assuming an initial thermal X-abundance (Ni = 1) 
and vanishing pre-existing asymmetry (N5_; = 0). The numerical results are compared 
with the analytic expression Eq. as indicated. 

The out-of-equilibrium decay regime is an efficient way to produce an asymmetry from 
the decays of heavy particles. However, it necessarily relies on the assumption that an 
initial X abundance was produced, either thermally or non-thermally, by some process at 
t => Tx and that one can neglect a possible Nj,_, generated during or after inflation and 


prior to the the onset of X-decays. Therefore, it is evident that such a scenario is plagued 
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by a strong dependence on the initial conditions and hence it requires to be complemented 
with a picture able to specify them, for example a detailed description of the inflationary 


stage. 


2.5 Inverse decays 


The case of out-of-equilibrium decays is valid rigorously only in the limit Kx — 0. If one 
defines zg as that special value of z corresponding to the time when the age of the universe 
is equal to Tx, i.e., t(za) = Tx, then, for Ky « 1, one has indeed zg ~ \/2/Kx de 
On the other hand, for Ky = 1, the bulk of X particles will decay at a temperature 


~N 


Ta = Mx/z2a = Mx/V2 so that inverse decays cannot be neglected as we did so far. The 


kinetic equations and are then generalized in the following way [?,}65} {63} 


X= —Dx(2) [Nx(2) - NBC) a 
a aa O 29) 


In the first equation, for Nx, the second term accounts for the inverse decays that, 
importantly, can now produce the X’s. On the other hand, one can see that a new 
term appears in the second equation for the asymmetry too, a wash-out term that tends 
to destroy what is generated from decays. This term is controlled by the (inverse decays) 
wash-out factor given by 


Nea 
Wee pe Sek 30 
X 2 X Ne X X, ( ) 


where n denotes the number of baryons or leptons in the X decay final state (as we will 
see n = 1 in the case of leptogenesis) and Noi is the equilibrium abundance either of 
leptons or baryons, in any case massless, produced from X decays. Note that the decay 
parameter K x is still the only parameter in the equations and thus the solutions will still 
depend only on Ky. They can be again worked out in an integral form [65]. In the case 
of the B — L asymmetry one can write the final asymmetry as 


= ie dz' WIP ( 


Np-1 = Np-re *) + ex Ke, (31) 


where now the efficiency factor is given by the integral 


oe dN. co Wt 1i 
k(Kx, 3) = -f dz’ pa ede ewe"), (32) 


13The equations and (29) are actually not only accounting for decays and inverse decays but also 
for the real intermediate state contribution from 2 + 2 scattering processes. This term exactly cancels 
a CP non conserving term from inverse decays that would otherwise lead to an un-physical asymmetry 
generation in thermal equilibrium (69). 
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In the limit Ky — 0 one has Wi}? — 0 in a way that the out-of-equilibrium case is 
recovered. In general, one can see that the wash-out has the positive effect to damp a 
pre-existing asymmetry but also the negative one to damp the same asymmetry generated 
from decays, thus reducing the efficiency of the mechanism. A quantitative analysis is 
crucial and it is very useful to discuss separately the regime of strong wash out for Ky = 1 


and the regime of weak wash-out for Ky < 1. 


2.5.1 Strong wash-out regime 


The strong wash-out regime is characterized by the existence, for Kx = 3, of an interval 
[Zin Zout] Such that WP 


Practically all the asymmetry produced at z < Zout is washed-out including, noteworthy, a 


= 1 and thus such that inverse decays are in equilibrium. 


pre-existing one. Moreover the calculation of the residual asymmetry is made very simple 
by the possibility to use the close equilibrium approximation, i.e., 
dNx daN% 2 


dz! dz “Kez Wx (2). eo) 


In this way the integral in Eq. can be easily evaluated for z; < 1 (64]|65). Indeed, 
this can be put in the form of a Laplace integral, 


he = f dz! et) | (34) 
0 


receiving a dominant contribution only from a small interval centered around a special 
value zp such that dw/dz = 0. In this way one can use the approximation of replacing 
WP (z") with WP (z") zg/2" in Eq. With this approximation the integral can be 
easily solved, obtaining 


2 n z 
kp a (1-7 a *). (35) 


For large Kx >> 1 and n = 2 this expression coincides with that one in The 
calculation of zg proceeds from its definition, (d/dz).,, = 0, approximately equivalent 


to the equation 


WP(n) = ey Ce (36) 


ZB 


14Tt is analogous but not coinciding with a saddle point approximation since in this case one does not 


Taylor expands the function 7)(z’, z) about the minimum. 
Note, however, that the definition for Ky has to be used instead of Ky = (1/2) (T p/H)z=1. 
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Figure 2: Final efficiency factor as a function of the decay parameter Ky for thermal 
and dynamical initial Nx abundance. The solid lines are the numerical solutions of the 
Eq.’s and (29), the short-dashed lines are the analytic results (see Eq. and 
Eqs. (44)+(45)), the dot-dashed line is the power law fit Eq. (141). The dashed box is 
the range Ky = 8-45 favoured by neutrino mixing data in the case of leptogenesis. 


This is a transcendental algebraic equation and thus one cannot find an exact analytic 


solution (see for an approximate procedure). However, the expression 
x 0.13 o- ni 
zg(K) ~2+4(n Ky)" e "5x, (37) 


provides quite a good fit that can be plugged into Eq. thus getting an analytic 
expression for the efficiency factor. For very large Ky this behaves as a power law kẹ x 
Ky’. In Fig. 2 we show a comparison of the analytic solution for «sẹ Eq. with the 
numerical solution (for n = 1). One can see how for Kx = 4 the agreement is quite good 
(at the percent level). Note that the Eq. implies that for large values of Ky one 
has zg & Zout, that particular value of z corresponding to the last moment when inverse 


decays are in equilibrium (WIP? > 1). In this way almost all the asymmetry produced for 
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Figure 3: Comparison between analytical (short-dashed lines) and numerical (solid lines) results 
in the case of strong wash-out (Kx = 100) for |ex| = 0.75 x 107° (from (64]). 


z S zp is washed-out and most of the surviving asymmetry is produced in an interval just 
centred around Zout, simply because the X abundance gets rapidly Boltzmann suppressed 
for z > Zou. An example of this picture is illustrated in Fig. 3 for Kx = 100 (from (64}). 
Instead of the abundance Nx, we plotted the deviation from the equilibrium value, the 
quantity A = Ny — NS’. The deviation grows until the X’s decay at z œ zg, when it 
reaches a maximum, and decrease afterwards when the abundance stays close to thermal 
equilibrium. Correspondingly, the asymmetry grows for z < 2a, reaching a maximum 
around z œ 1, and then it is washed-out until it freezes at zg œ Zou. The evolution of 
the asymmetry Ng_,(z) can induce the wrong impression that the residual asymmetry is 
some fraction of what was generated at z ~ 1 and that one cannot relax the assumption 
zi < 1 without reducing considerably the final value of the asymmetry. However, what is 


< zp there is no much 


~N 


produced for z S Zout is also very quickly destroyed and as fas as 2; 
dependence of the final asymmetry on z;. A plot of the quantity w(z, 00), as defined in 
the Eq. and shown in Fig 4 (from (64)), enlightens some interesting aspects. This is 
the final asymmetry that was produced in a infinitesimal interval around z. It is evident 


how just the asymmetry that was produced around zg survives and, for this reason, the 
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Figure 4: The function w(z’, co) for Ky = 10 (solid line) and Kx = 100 (dashed line, from lea). 


temperature Tg = Mx/zg can be rightly identified as the temperature of baryogenesis for 
these models. It also means that in the strong wash out regime the final asymmetry was 
produced when the X particles were fully non-relativistic implying that the simple kinetic 
equations and (29), employing the Boltzmann approximation, give actually accurate 
results and many different types of corrections, mainly from a rigorous quantum kinetic 
description including thermal effects, can be safely neglected. 

This is not the only nice feature of the strong wash-out regime. Since any asymmetry 
generated for z S zg gets efficiently washed-out, one can also rightly neglect any pre-existing 
initial asymmetry Nj,_,. At the same time the final asymmetry does not depend on the 
initial X abundance. In Fig. 5 we show how even starting from a vanishing X abundance, 
the X’s are rapidly produced by inverse decays in a way that well before zg the number of 
decaying neutrinos is always equal to the thermal number. As we said, the final asymmetry 
does not even depend on the initial temperature as far as this is higher than ~ Tg and 
thus if one relaxes the assumption z; << 1 to zi S zg — Azz, the final efficiency factor gets 
just slightly reduced (for example for Azg ~ 2 this is reduced approximately by 10%). 

Summarising, we can say that that in the strong wash out regime the reduced efficiency 
is compensated by the remarkable fact that, for Ti = Tg, the final asymmetry does not 


depend on the initial conditions and all non relativistic approximations work very well. 
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Figure 5: Fast thermalization of the X abundance in the strong wash-out regime. The final Ng_ rz 
abundance (for ex = 0.75 x 1076) is the same as in the case of an initial thermal abundance (cf. 


Fig. 3) and it is independent on the evolution at z « zp. 


These conclusions change quite drastically in the weak wash-out regime. 


2.5.2 Weak wash-out regime 


For Kx <1, zg rapidly tends to unity (see Eq. (37)). In Fig. 2 the analytic solution for 
the efficiency factor, (see Eq. (35)), is compared with the numerical solution. It can appear 
surprising that, in the case of an initial thermal abundance, the agreement is excellent 
not only at large Ky = 4, but also at small Ay < 0.4, with some appreciable deviation 
only in the range 0.4 S Ky X 4. The reason is that when the wash-out processes get 
frozen, the efficiency factor depends only on the initial number of neutrinos and not on 
its derivative and thus the approximation Eq. introduces a sensible error only in the 


transition regime Ky ~ 1. 
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Eq. can be easily generalized to any value of the initial abundance if one can 
neglect the X particles produced by inverse decays. More generally, one has to calculate 
such a contribution and it is convenient to consider the limit case of an initial vanishing 
X abundance. The production of the X’s lasts until z = 2.4, when its abundance is equal 


to the equilibrium value, such that 
Nx (2eq) = NX (Zea) - (38) 


At this time the number of decays equals the number of inverse decays. For z < Zeg, 
decays can be neglected and Eq. becomes 


dN. 
FE = Dx) NEC): (39) 
z 
For z < 1, one then simply finds 
K 
Nx(z) = = (2° — 28) . (40) 


In the weak wash-out regime the equilibrium is reached very late, when neutrinos are 
already non relativistic and Zeęą >> 1. In this way one can see that the number of X 


particles reaches, at z œ 2g, a maximum value given by 


Ny (za) = N(Kx) = Kx. (41) 
For z > Zeq inverse decays can be neglected and the X’s decay out of equilibrium in a 
way that 
Nyx(2z > zea) @ Nx(zeale Fea 2) DO (42) 
Let us now consider the evolution of the B—L asymmetry calculating the efficiency factor. 
Its value can be conveniently decomposed as the sum of two contributions, a negative one, 
kr , generated at z < Zq, and a positive one, k}, generated at z > Ze; In the limit of 
no wash-out we know that the final efficiency factor must vanish, since we have seen that 
in this case xp = Ni and we are assuming a vanishing initial abundance. This implies 
that the negative and the positive contributions cancel with each other. The effect of the 
wash-out is to suppress the negative contribution more than the positive one, in a way 
that the cancellation is only partial. In the weak wash-out regime it is possible in first 
approximation to neglect completely the wash-out at z > 2% . In this way it is easy to 
derive from the Eq. the following expression for the final efficiency factor: 


ke = N(Kx) — 2 (1 — emi neo) : (43) 


One can see how it vanishes at the first order in N(Kx) « Kx and only at the second 
order one gets Ky ~ (917/64) KX. 
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2.5.3 Final efficiency factor: global solution 


Generalising the procedure seen for the strong wash-out, it is possible to find a global 
solution for K¢(/Kx) valid for any value of Kx. The calculation proceeds separately for 


«& and «* and the final results are given by 
Ke (Kx) S=2 i NE) Cu - 1) (44) 


and 
2, 


ct(Ky) = 
nel x) zp(Ky)nKyx 


(1 = i) ; (45) 


The function N (Kx) extends, approximately, the definition of N(Kx) to any value of Kx 


and is given by 
N(Kx) 


Ji 
(1 +/N (Kx)) 
The sum of the Eq.’s and is plotted, for n = 1, in Fig. 2 (short-dashed line) and 


compared with the numerical solution (solid line). 


N(Kx) = (46) 


We can now draw a few conclusions about a comparison between the weak and the 
strong wash-out regime. A large efficiency in the weak wash-out regime relies on some 
unspecified mechanism that should have produced a large (thermal or non-thermal) X 
abundance before their decays. On the other hand, the decrease of the efficiency at large 
Ky in the strong wash-out regime is only (approximately) linear and not exponential 
[64]|65}. This means that for moderately large values of Kx a small loss in the efficiency 
would be compensated by a full thermal description such that the predicted asymmetry 
does not depend on the initial conditions, a nice situation that resembles closely the 
situation in standard big bang nucleosynthesis, where the calculation of the primordial 
nuclear abundances is also independent of the initial conditions. 

Another point to notice is that we discussed a very general picture taking into account 
just decays and inverse decays. However, considering specific models, one might have to 
consider some specific processes. We will discuss in Section 5 the case of leptogenesis 
and we will see that in that case scatterings also need to be taken into account for the 
production of decaying particles producing the asymmetry, though in the strong wash-out 
regime this will not have a great impact on the final predicted asymmetry. On the 
other hand, we will see that flavour effects can dramatically change the calculation of the 


asymmetry in certain situations. 
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2.6 GUT baryogenesis 
The advent of GUT theories in the seventies provided a very well motivated 


realistic extension of the SM to embed baryogenesis from heavy particle decays. Here we 
want briefly to show how the general results specialise considering a simple toy model} 

The asymmetry is in this case generated by the baryon number violating decays of 
super heavy bosons assumed to be complex scalars. Since we need interference of tree 
level and one-loop diagrams for the total CP asymmetries not to vanish, one needs at 
least two of them that we can denote them by X and Y and their masses by Mx and 
My respectively. If we assume My >> My and Ky > 1, we can neglect the asymmetry 
produced by the Y’s in a way that the B — L asymmetry will be dominantly produced 
by the X’s and their anti-particle states X. Let us assume that there are two baryon 
violating decay channels|!"| X + Yı +yz and X > Y, + Ya, and its CP conjugated, 
X > yı +yz and X > 3 + 4, where each y; (i = 1,...,4) is a massive fermion with 
baryon number B;. 

In this case the kinetic equations and still hold but now Nx should refer to 
the total abundance of X and X particles and the washout factor is given by Eq. (30) 
with n = 2. 

In realistic models the masses of super heavy gauge bosons is close to the grand-unified 
scale Maur. Since the proton lifetime is approximately given by Tp ~ Onur M$ my >, the 
current experimental lower bound 7, = 10** years translates into very stringent lower 
bounds, Mx = 10!°-'° GeV, for the masses of the superheavy gauge bosons that imply 
correspondingly very high reheat temperatures in the early universe in tension with upper 
bound from CMB anisotropies. We will see that SO(10) GUT theories also predict the 
existence of RH neutrinos that can be much lighter. Their decays can then be responsible 
for the production of the B — L asymmetry and this is at the basis of leptogenesis that 
we will discuss in detail in Section 5. Therefore, from this point of view, leptogenesis 
can be regarded as a way to achieve successful baryogenesis at lower reheat temperatures 
within GUT models with the addition of RH neutrinos. Moreover, the introduction of RH 
neutrinos also provides an attractive solution to the explanation of light neutrino masses 
and mixing and, therefore, leptogenesis somehow appears as the most likely way how the 


asymmetry can be produced within these models in a thermal way. 


16For a more complete discussion we refer the reader to |65}. 
17The existence of at least two channels is crucial since X is a boson, and in this case the CP conjugated 


state would decay exactly with the same rate because of CPT invariance. In the case of leptogenesis we 
will see that the decaying particle is a Majorana fermion and in that case one decay channel is sufficient, 


from this point of view it is more economical. 
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However, one can also consider a non-thermal production of the super heavy gauge 
bosons and in this case one can successfully reproduce the asymmetry for even much 
lower values of the reheat temperature. Various proposals have been made, such as 
production from inflaton decays [71], at preheating or from primordial black hole 
evaporation (24]/37|/38]/76}{82). 

Non-thermal production mechanisms can of course be considered also in general, 
beyond GUT baryogenesis, i.e., for some generic heavy particles producing the asymmetry. 
In this respect recently it has been also considered the case of non-thermal heavy particle 


production from collisions of runaway bubbles originated during first order phase transitions 


384) 


3 Dark matter 


Very early hints of the existence of a dark matter component from anomalous motion of 
stars in our galaxy were already discussed at the beginning of the twentieth century and 
the same name, matiére obscure in French, was already introduced by Poincaré in 1906 
[85]. On much large scales, the anomalous dynamics of cluster of galaxies in 1930’s 
also pointed to the existence of such component. A more solid evidence was found only 
much later in the 1970’s in galactic rotation curves [87][88]. An understanding of the large 
scale structure of the universe then provided the first evidence of the non-baryonic nature 
of the dominant component of dark matter and of its primordial origin, that had to occur 
prior to the matter-radiation equality time teq ~ 52, 000 yr. 

It also became clear at the beginning of the 1980’s that such component could not be 
constituted by massive neutrinos, since still relativistic and free-streaming too fast at teq 
(hot dark matter), but rather by some new kind of non-standard matter that had to be 
already non-relativistic at teq (cold dark matter). It was then only with the discovery of 
CMB acoustic peaks that it was possible to firmly establish its non-baryonic nature that 
cannot be explained within the SM. It is to this non-standard component that one usually 
refers simply as dark matter. The cold dark matter contribution to the energy density 
parameter is today well determined both from dynamics of clusters of galaxies and from 
CMB anisotropies and from Planck + BAO data it is found 


OQpmoh” = 0.11933 + 0.00091 (47) 
and from this, using Eq. for the value of h, one finds 


Qpmo = 0.261 + 0.005. (48) 
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The total matter contribution, baryonic plus cold dark matter, is then given by 
Ooh? = 0.14240 + 0.00087 (49) 
and using again Eq. for h one finds 
Qmo = 0.311 + 0.006. (50) 


The most popular idea is that dark matter is made of new (elementary or composite) 
particles whose existence cannot be explained within the SM|® Currently, the discovery of 
gravitational waves from black hole merging has drawn a lot of attention on the alternative 
option that dark matter is made of primordial black holes with masses that can be 
potentially as large as hundred solar masses, though a host of observational constraints 
rules out the possibility of a population of monochromatic primordial black holes making 
the whole dark matter, except for a window 3.5 x 1071" S Mppu/Mo Ș 4 x 10°” 
corresponding to asteroid masses. The possibility of a population of primordial black 
holes with a spread mass spectrum is, however, not yet completely excluded p3]. Another 
possibility is that dark matter can be actually understood in terms of models of modified 
gravity rather than as an additional matter component, the most popular proposal is 
the MoND theory or its covariant TeVeS theory [94]. Recently mimetic dark matter has 
also attracted attention since it also offers a unified solution to the dark energy puzzle 
and even a model for inflation [95]. However, thorough statistical analyses including all 
cosmological observations disfavour such alternative models compared to the standard 
ACDM model or variations to it where dark matter is assumed to be made of particles. 


We will then focus our attention on particle dark matter models. 


3.1 Particle dark matter models 


In particle dark matter models, in addition to the requirement that the final relic abundance 
reproduces the measured dark matter contribution to the total cosmological energy density 
in Eq. (47), the dark matter particles need to be long-lived on cosmological scales. In 
order to reproduce galactic rotation curves and the dynamics of clusters of galaxies, the 
life-time of dark matter particles Tp), has to be necessarily longer than the age of the 


Universe to, so that one has to impose a lower bound 


Tom Z to 4x 10'"s. (51) 


18Recently the possibility that dark matter is made of (SM) uuddss hexa-quarks has been explored 
in [89], showing that a mass ~ 1.2 GeV would be required but that on the other hand this is ruled out 
by the stability of Oxygen nuclei. 
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The possibility of dark matter particles decaying with Tom œ to has been seriously 
considered to solve well known problems of the standard CDM scenarios in reproducing 
observations on scales equivalent to that of dwarf galaxies such as the missing satellite 
problem and the too big to fail problem. This decaying dark matter scenario is however 
currently disfavoured compared for example to other scenarios that aims at solving the 
same problems, in particular models with self-interacting dark matter. It should also 
be said that in recent years with more accurate astronomical observations and N-body 
simulations that include baryon feedback, the shortcomings of the CDM scenario tend to 
be greatly ameliorated and currently there is no compelling motivation to abandon it [| 

If the dark matter particle decays into SM particles, then the lower bound becomes 
many orders of magnitude more stringent and the exact value depends specifically on the 
mass of the dark matter particle and on its decay products. We will discuss specific cases. 
More generally, the coupling of dark matter particles to SM particles has to be weak 
enough to escape all experimental constraints and of course also that DM particles need 
to decouple prior to the onset of big bang nucleosynthesis. For this reason dark matter 
particles cannot interact electromagnetically and strongly (though their constituents can 
in case of composite dark matter particles). 

The dark matter number density of dark matter particles at present is related to the 


energy density parameter simply by 


Npmo h? 
A OOM ajg t A, (52) 


Mpm DM 


NpMo = Eco h` 


From this expression we can also derive an expression for the abundance of dark matter, 
relatively to that of photons, at the production and more precisely at the freezing time 


tp. Assuming that entropy is conserved between freezing and present time, one finds 


N pumo h? 
( 2 n a O a (53) 
N; J: Nyo gso no Mom 9so Mom 


where in the last numerical expression we assumed that T; = T (tt) Z 100 GeV so that we 
could use the SM value g2“ = 427/4 for the number of ultra-relativistic degrees of freedom. 
This simple expression gives some insight on the efficiency of the production mechanism. 
In particular, it shows that for higher masses we need a less efficient mechanism since one 


needs a smaller number density. If entropy is not conserved, then one has to include a 


19There are of course other tensions in the ACDM model, the strongest being the Hubble tension, but 
the decaying dark matter scenario still does not help addressing them. Solutions with some modification 
of the ACDM model in the pre-recombination stage seem to be the most favoured from this point of view 
(for an extensive discussion see (96)). 
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further factor So/S;: > 1, showing that since the relic abundance gets additionally diluted 
compared to photons by an entropy increase, one needs a more efficient mechanism of 


dark matter production to compensate. 


3.2 The WIMP paradigm 


We mentioned that massive ordinary neutrinos play the role of dark matter but of the 
unwanted (hot dark matter) type. For this reason, one gets a stringent upper bound on 
the sum of the neutrino masses. However, if one considers the case of much heavier neutral 
particles with weak interactions, then a new solution appears. In the case of light active 
neutrinos they are produced thermally and they decouple at a temperature TY. ~ 1 MeV 
when they are still ultrarelativistic. If one considers considers a weakly interacting massive 
particle (WIMP) X, its mass is assumed in the range mx ~ 10GeV-1 TeV, since it is 
typically associated to models where the scale of new physics is just above the electroweak 
scale to address the naturalness problem. In this way, having weak interactions, they 
decouple when they are fully non-relativistic at a freeze-out temperature T ~ mx /xs, 
with x ~ 25. This can be easily estimated using a simple instantaneous decoupling 
approximation, imposing the usual criterion .y(T;*) ~ H(T;*‘) 97|. 

A more accurate result can be derived solving a rate equation”) for the number density 
of WIMPs nx(t), the Lee-Weinberg equation 

dnx 


ae +3 Hnx(t) = — (Tann Bret) [nx (t) E (nD H) ) (54) 


where (Cann Pret) is the thermally averaged annihilation cross section times the Moller 
velocity. Solving analytically this equation, one finds that the final contribution to the 


energy density parameter from WIMPs is 


4x 107 (x2. 1.6 x 10° om? | 
(Tann Brel) £ (Tann Breit 


where (dann Brel) is calculated at the freeze-out time. Typical weak values of thermally 


xoh? œ= 


(55) 


averaged cross sections are given (order-of-magnitude-wise) by 


2 100 GeV \ ? 
loaba aea E a aa { (56) 
f 2 
my mx 
2 
~ 4x 107?" cm? 100 GeV 
Mx ` 


20Tt can be derived from the Boltzmann equation for the X particle distribution function integrating 


over momenta and within certain approximations such as the Maxwell-Boltzmann approximation. 
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Here we used aw ~ 0.03, the value of the dimensionless weak coupling constant at energies 
~ 100 GeV and for the relative velocity at the freeze-out a typical value p,e œ 0.1. It is 
then intriguing that for WIMP masses mx ~ 100 GeV-—1 TeV, the expected values from 
naturalness, one obtains Qyxoh? ~ 0.1, nicely reproducing the observed value of the dark 
matter abundance. Because of this tantalising WIMP miracle, for long time WIMPs have 
been regarded as the most attractive dark matter candidate also in consideration of their 
detectability by virtue of their weak interactions. 

Some of the proposed most popular realistic examples of WIMPs |100!101], are associated 
with specific solutions to the naturalness problem. Supersymmetry is certainly the most 
popular solution, perhaps also the most elegant one, and for this reason the most popular 
dark matter WIMP candidate is the lightest supersymmetric neutral spin 1/2 fermion, 
the so-called neutralino. This can be made stable by adding a discrete symmetry called 
R-parity with an associate conserved quantum number with only two discrete values: +1 
(even particles) or -1 (odd particles). Supersymmetric and standard model particles have 
opposite R-parity, minus and plus, respectively. Therefore, a supersymmetric particle 
cannot decay just into SM particles and this implies that the lightest supersymmetric 
particle has to be stable. Neutralino as dark matter WIMP was proposed in the eighties 
and it became such an attractive solution that its discovery was considered 
just a matter of time. Moreover, as we discussed in Section 2.3, within supersymmetry 
one can also have successful electroweak baryogenesis, thus realising a unified picture for 
the origin of matter in the universe "| Unfortunately, intense searches of different kinds 
have not fulfilled expectations and this has stimulated new ideas for models solving the 
naturalness problem with their own associated dark matter WIMP candidate. 

A popular class of theories providing a solution to the naturalness problem, alternative 
to supersymmetry, are theories with extra-dimensions. They also usually contain candidates 
for (thermally produced) dark matter WIMPs, the most popular one being the lightest 
Kaluza—Klein particle. This is associated to some standard model particle, typically the 
hypercharged gauge boson B. Other recent examples of models addressing the naturalness 


problem are (in brackets their associated dark matter WIMP): 
(i) little Higgs models (T-odd particles); 
(ii) two Higgs doublet models (neutral Higgs boson); 


(iii) twin Higgs models (twin neutral Higgs); 


2lThis possibility has been investigated in particular within the MSSM in the presence of a light 


stop and within the nMSSM [t06]. 
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(iv) dark sector models (dark lightest particle). 


The evidence found so far for the existence of dark matter is uniquely based only 
on its gravitational effects. The attractive feature of dark matter WIMPs is that they 
also interact weakly and this implies specific experimental signatures that allow to test 
the WIMP paradigm. There are three strategies pursued to detect dark matter WIMPs 
through their weak interactions [101]/107}. If we indicate with SM some generic standard 
model particle, and with X our WIMP, one looks for signals of the following kind of 


interactions: 


(i) X+SM—>X+SM (direct searches) ; 


(ii) X+X —> SM +5SM (indirect searches) ; 
(iii) SM + SM > X’s + SM’s (collider searches) 


Direct searches. Dark matter WIMPs, though very rarely, can collide with ordinary 
matter through their weak interactions. In particular, they can hit nuclei in a detector 
and transfer them energy (nuclear recoil energy) that can be detected. In the last three 
decades different experiments have been placing more and more stringent limits on the 
mass and on the cross section with nuclei of the dark matter WIMPs. However, there are 
also models of WIMPs with spin dependent cross section and in this case constraints are 
weaker. These constraints also depend on the evaluation of astrophysical quantities such 
as the local density and the velocity distribution of WIMPs at our position in the Milky 
Way, clearly affected by some theoretical uncertainties. 

Indirect searches. A complementary strategy is to look for the products of dark matter 
interactions within astronomical environments. The interactions with ordinary matter are 
in this case unhelpful. In a perfectly homogeneous universe, dark matter annihilations, 
though not completely turned off, would be so strongly suppressed that it would be 
impossible to detect any observable effect. However, dark matter, like ordinary matter, 
clumps on scales smaller than 100 Mpc. Since the rate of annihilations x pôw, one 
expects that the flux of annihilation products is greatly enhanced in very dense dark 
matter regions, such as the galactic centre or dwarf galaxies, where the density of dark 
matter particles is expected to be order of magnitudes above the average value. First of 
all one can hope to detect y-rays coming from such dense dark matter regions. These 
have the advantage that would travel from the production site to us in a straight line, 
retaining information on the source position in the sky. Therefore, one can target dense 


dark matter regions to have better chances to distinguish a signal from the background. 
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Ideally, one would like to observe a monochromatic line but unfortunately the rate of 
WIMPs direct annihilations into photons is too weak and what one can realistically 
detect are photons produced as secondary particles from the primary particles produced 
in the annihilations. These generate an excess spread on a wide range of energies. In 
the mass range (0.1-1) TeV the most stringent upper bound on the WIMP annihilation 
cross section has been placed by the Fermi-LAT space telescope and it basically excludes 
the typical values expected for thermally produced WIMPs the value at freeze-out in 
Eq. (56). At values of the mass above TeV, most stringent constraints are placed by 
the HESS ground based y-ray telescope but they still allow values (dann ()¢ as large as 
10735 cm?. On the other hand at values my < 0.1 TeV the most stringent limits come from 
the Planck satellite observations of CMB anisotropies constraining (Cann Brei) at the time 
of recombination. These constraints strongly exclude weak cross sections for thermally 
produced WIMPs. 

Alternatively, one can search for dark matter products in charged cosmic rays. However, 
in this case particles are deflected by the galactic magnetic fields and would travel to us 
along a random curvy path inside the galaxy, so that directionality is lost. One can 
simply compare the measured spectrum to the one predicted by the SM looking for some 
excess that would be the result of WIMPs annihilations. Typically, instruments detect 
electrons and protons, that can originate from many astrophysical sources, and above 
all positrons and anti-protons [109], that is even more interesting for dark matter 
searches since these would be produced in equal amounts to electrons and protons in dark 
matter annihilations. Given a source for cosmic rays with a specified energy spectrum 
at production, the prompt spectrum of particles, the predicted energy spectrum on the 
Earth is the result of many complicated processes, with the main one being the diffusion 
due to the galactic magnetic fields. This is usually calculated solving a Fokker—Planck 
equation [110], that is a generalised diffusion equation. Cosmic rays are detected with 
balloon-type, ground based telescope arrays or satellite-based detectors. In 2009 the 
PAMELA satellite detector reported an excess in the positron spectrum in the energy 
range (20-200) GeV. The Fermi-LAT and the AMS-02 satellite detectors confirmed the 
excess though to a lower level. A dark matter interpretation seems plausible but it 
encounters great difficulties since the value of the cross section required to explain the 
excess is much larger than the typical values expected for thermally produced WIMPs 
(see Eq. (56)). Many specific models have been proposed but they are in tension with the 
constraints from y-rays so that further experiments are needed for a conclusive verdict. 

Neutrinos can also be very useful for indirect searches of WIMPs. Thanks to their 


weak interactions, WIMPs would dissipate energy scattering off nuclei in central dense 
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regions of celestial bodies such as stars and planets. If their velocity becomes smaller than 
the escape velocity, they get trapped accumulating gradually in the central dense regions 
with a capture rate Deap¢. When density gets higher and higher, the annihilation rate Dann 
becomes sufficiently large to balance the capture rate and equilibrium is reached. WIMPs, 
such as neutralinos, cannot annihilate directly into neutrinos but their annihilation products, 
e.g., quarks and gauge bosons, would produce them secondarily. Neutrinos are the only 
particles able to escape the centre of stars travelling in a straight line and retaining 
information on their production site. In our case the only object able to produce a 
detectable neutrino flux would be the Sun or, to a lower level, the same centre of the 
Earth. Neutrino telescopes, big neutrino detectors able to track the arrival direction 
of neutrinos, pointing toward the centre of the Sun or the Earth should then detect a 
high energy neutrino flux. This strategy is sensitive both to the scattering off nuclei and 
annihilation cross section of WIMPs. The derived upper bounds on spin independent 
and on annihilation cross section are less stringent than those derived from direct and 
y-rays searches. However, those on spin dependent cross section, derived by combining 
data from ANTARES, IceCube and SuperKamiokande neutrino telescopes, are the most 
stringent ones. 

Collider searches. Dark matter WIMPs could also be produced in colliders and 
currently the LHC has the right energy to test a wide range of masses. Unfortunately, so 
far, there are no signs of new particles in the LHC and not even of missing energy and 
momentum that cannot be explained with neutrinos. From LHC results one can then 
place upper bounds on the dark matter production cross sections and these are perfectly 
compatible with those from direct and indirect searches. 

We can fairly conclude that current experimental results rule out simple ‘WIMP 
miracle’ expectations, strongly motivating modifications or alternative models. 

The main reason why the WIMP miracle is ruled out is that the upper bounds on 
the cross sections are so stringent to point to much smaller values than those required 
by Eq. (56): from Eq. one would then obtain a dark matter abundance much higher 
than the observed one. However, one can still think to save the idea of WIMP dark 
matter relaxing one or more assumptions on which the WIMP miracle scenario relies 
on. For example, one assumption is that the freeze-out occurs in isolation, i.e., that it 
is not influenced by the existence of possible additional new particle species that could 
also be WIMPs. Relaxing this assumption, one can obtain the correct abundance of 
dark matter WIMPs, with much lower values of (o 3)s, taking into account three possible 
effects [111]. A first important one is given by so-called co-annihilations, occurring when 


the mass of the dark matter WIMP is sufficiently close to the mass of a heavier unstable 
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WIMP. Another possibility is that annihilations occur resonantly, at energies close to the 
Z or Higgs boson mass, and this also tends to enhance the annihilation rate reducing 
the final relic abundance compared to the non-resonant case and making it compatible 
with the observed value. Both effects can be realised within supersymmetric models, 
singling out particular regions in the space of parameters (the so-called co-annihilations 
and funnel regions) where the neutralino is the lightest supersymmetric particle and its 


relic abundance reproduces the observed dark matter abundance. 


3.3 Beyond the WIMP paradigm 


The list of models beyond those of dark matter WIMPs relying on a traditional freeze-out 
production mechanism is impressively long. This certainly shows how the dark matter 
puzzle is one of the greatest challenges in modern science. Here we mention the most 
popular ideas referring the reader to more specialistic reviews (101/112). 


e Hidden dark matter and the WIMPless miracle. The WIMP miracle can be somehow 


revisited considering dark matter particles with interactions such that 


a? MWIMP WIMP 
(Tann Brel) t = m a (Tann Bye ; (57) 
X WwW 


where (Cann 8) PP denotes the special value of (Cann Breie, given by Eq. (56), for 
mx = mwmpr. For ax/mx = aw/mwimp one still has (ann Breit = (Cann Bret)”. 
This means that even though the dark matter particle is not a WIMP, it has the 
same annihilation cross section of WIMPs. In this way the observed dark matter 
abundance is still reproduced through the usual thermal freeze-out mechanism. For 
example, a dark matter particle species with mx ~ MeV and a ~ 10-° aw would 
still yield the correct relic abundance but it would escape the tight experimental 


constraints we discussed. 


e Feeble, extremely or super-interacting massive particles (FIMPS or SWIMPs) 
114|. These dark matter particle candidates have interactions much weaker than 
weak interactions. Nevertheless, they can be still produced thermally. However, 
their abundance never reaches the thermal value before freeze-out, rather it directly 
freezes-in to the correct value starting from an initial vanishing abundance [115]. 
The relic abundance is in this case proportional to the annihilation cross section 
instead of inversely proportional. In this way one can obtain the correct abundance 


for the same values of masses as in the case of WIMPs but with much smaller cross 
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sections, thus escaping current experimental constraints. Sterile neutrinos can be 
regarded as a particular important example of FIMP dark matter, in fact the first 


to be proposed. We will discuss specific models both in Section 6 and Section 7. 


Alternatively, their production can occur non-thermally from late decays of heavier 
particles. For example, one could have an unstable WIMP that, by virtue of the 
WIMP miracle has the correct relic abundance Qwmyp. Even though these unstable 
WIMPs cannot play the role of dark matter, still their abundance can be transferred 
to a Super-WIMP (cosmologically stable) particle through decays, in such a way that 
OQswime = Qwimp os (58) 
TWIMP 
If mswmip is not too much smaller than mwmp, the WIMP miracle still works 
also for SWIMPs, with the difference that SWIMPs escape direct and indirect 
search constraints. However, WIMPs could be still produced if some astrophysical 
environments and their decays have some testable effects in cosmic rays if the 
decaying WIMP is charged. Particle colliders may also find evidence for the SWIMP 
scenario in this case. Finally decays of WIMPs into SWIMPs might have an impact 
on CMB and BBN in the form for example of extra-radiation, sometimes called 
dark radiation and often parameterised in terms of the number of effective neutrino 
species. In the case of BBN, for life-times between ~ 1s and ~ 1000s, the products 
of decays alter in general the values of the primordial nuclear abundances in an 
unacceptable way yielding constraints on the parameters of the model. However, 
for some fine-tuned choice of the parameters, the late WIMP decays could even be 
beneficial, solving the long-standing lithium problem in BBN [116]. 


A popular example of SWIMP particle is the gravitino, superpartner of the graviton 
with spin 3/2. It can play the role of dark matter when it corresponds to the lightest 
supersymmetric particle, since in this case it can be cosmologically stable. It can 
be produced both thermally, through the freeze-in mechanism, and non-thermally 
through the decays of the next-to-lightest supersymmetric particle (typically the 
neutralino). In order to get the correct abundance, the reheat temperature cannot 
be above ~ 10!°GeV. For gravitino masses below ~ 1TeV, this upper bound on 
the reheat temperature, becomes a few orders of magnitude more stringent taking 
into account the constraints from BBN. Such stringent upper bound on the reheat 
temperature is incompatible with lower bounds on the reheat temperature within 
different scenarios, thermal GUT baryogenesis that we discussed and, as we will see, 


even with the more relaxed analogous lower bound arising in thermal leptogenesis: 
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this is the well known gravitino problem [117]. 


Azions. Another attractive and very popular dark matter candidate emerges naturally 
from the so-called Peccei-Quinn mechanism for the solution of the strong CP violation 
problem in QCD [118]. This mechanism predicts the existence of a new particle, 
called azion, that can have the right features to play the role of a dark matter 
particle. Axion particles would be associated to a pseudoscalar field a with coupling 
to gluons determined by a coupling constant f7t, where fa has the dimension of 
energy, determining the decay rate of axions and also the same axion mass given 
approximately by 


(59) 


Ma ~ 107° eV (a ax) 


fa 
In order for axions to be stable on cosmological scales, one needs f, = 10° GeV but 
actually, from astrophysical constraints, one obtains a much more stringent lower 
bound, fa = 10° GeV translating into m, 


of axions would proceed non-thermally through a vacuum misalignment mechanism 


< 10meV. The cosmological production 


~~ 


producing an abundance 


Qa ~ 0.36? a = (60) 
: 102 GeV f 


where @ is the initial vacuum misalignment angle. The correct dark matter density 
is then achieved for ma ~ 107? 6° meV ~ (1076-1074) eV. Despite the small mass, 
axions would behave similarly to cold dark matter rather than hot dark matter, in 
agreement with structure formation constraints. This is because dark matter axions 


would basically form a Bose-Einstein condensate with negligible kinetic energy. 


The axion can be detected directly since it interacts with SM particles. Currently, 
the most stringent constraints on the mass and coupling constant are placed by the 
conversion of dark matter axions to photons through scatterings off a background 
magnetic field (Primakoff process). Within supersymmetry, the supersymmetric 
partner of the axion, the avino, is also a viable candidate for dark matter and, 


contrarily to the axion, it would behave as a SWIMP, similarly to the gravitino. 


Super massive dark matter (WIMPzillas). Many theories such as grand-unified 
theories, as we have seen when we discussed GUT baryogenesis, predict the existence 
of super massive particles with masses as high as the grand-unified scale, Mx ~ 
(1015-1016) GeV. Since they are above the maximum allowed value of the reheat 
temperature, Try < 10!° GeV, they cannot be produced thermally. However, they 
can be produced non-thermally at the end of inflation in different ways |119]: 
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(i) At reheating, when the inflaton field energy is transferred to all other particles. 
In this case the correct dark matter abundance can be attained even for masses 


that are orders of magnitude greater than the reheat temperature. 


(ii) By inflaton decays if the mass of the inflaton is greater than the mass of the 


supermassive dark matter particle. 


(iii) Finally, they can be produced at the so-called preheating, when the energy of 
the inflaton field can create supermassive particles through non-perturbative 
quantum effects on a curved space leading to parametric resonant production, 
an effect that can be mathematically described using Bogoliubov transformation 
leading to a Mathieu equation for particle production. Such a preheating 
stage at the end of inflation would also be responsible for the production of 


gravitational waves during inflation and in this way it could even be tested. 


3.4 Asymmetric dark matter 


Asymmetric dark matter modeldg”?] provide a very elegant solution to circumvent the 
experimental stringent constraints we mentioned for WIMPs but they also provide an 
attractive way to combine baryogenesis with dark matter models. Moreover they have 
the additional ambitious motivation to address the apparent coincidence Q mo ~ 5 Q po. 
The freeze-out mechanism assumes a vanishing initial asymmetry between the number 
of dark matter WIMPs and anti-WIMPs. However, as in the case of ordinary matter, some 
mechanism might have generated an asymmetry in the dark matter sector prior to the 
freeze-out. If this asymmetry is sufficiently large, then the relic abundance is basically 
given by the same initial asymmetry. This idea is quite attractive since one can also 
naturally link the problem of generation of the baryon asymmetry to the problem of 
dark matter, obtaining a combined model on the origin of matter in the universe. The 
indirect detection rate from annihilation is of course strongly suppressed compared to 
the standard case (since basically one is left only with particles or anti-particles today) 
and in this way the stringent constraints holding in the standard freeze-out scenario are 
evaded. On the other hand, signals in direct detection searches are usually enhanced and 
many constraints have been placed excluding different regions in the parameter space. 
In particular the elegant solution of having the same asymmetry in baryonic matter and 
dark matter, |n| = npm], implying mx = (Qpw/Qz)m, ~ 5 GeV, where m, ~ 1 GeV 


is the proton mass, is now ruled out by the stringent LUX constraints. 


22We just briefly discuss them here, we refer the reader to two recent reviews for an exhaustive discussion 


eae 
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4 Neutrino masses and mixing 


An explanation of neutrino masses and mixing necessarily requires an extension of the 
SM. It is then reasonable to search for scenarios explaining the origin of matter in the 
Universe that can be realised within those extensions of the SM also able to describe 
neutrino masses. In this section we review the current experimental status on neutrino 
parameters. In the next sections we will see how extensions of the SM that explain 


neutrino masses and mixing can also address the origin of matter in the universe. 


4.1 Low energy neutrino experiments 


SM neutrinos interact at low energies, when the SU(2), x U(1)y electroweak symmetry 


is broken to U(1)em, via neutral currents 


bu = bD Vii Yu VLi , (61) 


/=1,...,Ny 


where neutrinos couple to the Z bosons, and via the leptonic charged current 


J = (IE = he we: (62) 


where neutrinos couple to the W* bosons. Like all SM lepton interactions with gauge 
bosons, they are flavour blind, thus respecting lepton universality. Therefore, in the 
previous expression leptons fields have to be meant in a generic flavour basis. Flavour 
dependence originates from charged lepton Yukawa interactions. In a generic flavour basis, 


denoted by primed Latin indexes, these are described by the Lagrangian term 
-Lý = Ly hpp lpg ® + h.c., (63) 


where © = (+, 4°)? is the Higgs doublet. After EW symmetry breaking by a non zero 
vacuum expectation value (VEV) v of the Higgs field, the Yukawa coupling matrix hf j 


yield the charged lepton mass matrix 
M51 =i hing (64) 
and the charged lepton Lagrangian mass term can be written as 
-LE = Try mb y bg (65) 


In the charged current in Eq. (62), the charged lepton fields are not in general mass 


eigenfields. However, one can always diagonalise the charged lepton mass matrix by a 
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bi-unitary transformation of the left- and RH charged lepton fields, 


Lig ->ar = lra = Ut a liw, (66) 


lri > AR = Ra = Ua Cri , (a = e, ft, T) (67) 


so that correspondingly 


Mir = (u; diag(me, My, M7) Un) (68) 

uy! 
Therefore, the matrices Uz and Up transform the charged lepton LH and RH fields 
respectively from the generic primed basis to the weak basis. After this transformation 


the leptonic charged current (see Eq. (62)) becomes 


‘an = Soar UF as Vj! - (69) 


at! 


In this way the basis where the the charged lepton mass matrix mf is diagonal identifies 
a privileged flavour basis, the weak interaction basis, providing a natural lepton basis 
where to describe physical processes, since it is automatically (and quickly) selected by 
kinematics. This basis is determined by a bi-unitary transformation of the left and RH 


charged lepton fields 


fry > ar = lro = UF ei! lty, (70) 
rit > AR = lRa = Ung Crit , (a = e, fl, T) (71) 
in a way that 
Mi = (uz diag(Me, Mp, M7) UR). ' (72) 
age 


In the weak interaction basis the leptonic charged current can be written as 


Jee = SO Ub aw Ui (73) 


a,i! 


Defining the neutrino weak interaction eigenfields as 
VLa = UF eit Vii! , (74) 
the leptonic charged current can then be written as 


at = or VLa ,; (75) 
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showing that the neutrino weak interaction eigenfields create the neutrino states produced 
in the decay of a W boson in association with a physically charged lepton a of definite 
mass. 

In the SM neutrinos are massless. If neutrinos do have a mass, as now established by 
neutrino oscillation experiments, then one can define a basis of N neutrino mass eigenfields 
Vi, Va ..., Vy, associated respectively to masses m1, M2,..., My, that in general do not 
coincide with the weak eigenfields. The leptonic mixing matrix U is defined as the matrix 
that operates the transformation from the mass eigenfields v; to the weak eigenfields va, 


explicitly 
N 
y= y u% Vi. (76) 
i=1 


In general the number N of mass eigenfields mixing with the weak eigenfields does not 
have to coincide with the number of weak eigenfields (three) but it could be higher. This 
can well happen since, in extensions of the SM incorporating neutrino masses, in addition 
to three active neutrino fields there can also be sterile neutrinos mixing with the ordinary 
ones. In this case if one describes the mixing accounting only for three mass eigenfields, 
the (3 x 3) leptonic mixing matrix U,; with i = 1,2,3, also referred to as the PMNS 
matrix, is in general non-unitary. We will discuss a specific case with N > 3 within the 
type-I seesaw extension of the SM and we will see that the active-sterile neutrino mixing 
can have important cosmological applications especially in connection to the origin of 
matter. For the time being let us restrict ourselves to the standard case with N = 3 and 
unitary U. 

We can also consider a unitary transformation that brings from the neutrino mass 


basis to a generic primed basis: 
7 = Ufi VLi. (77) 


In this generic primed basis the leptonic mixing matrix can be regarded as the combination 
of two unitary transformations: a transformation U% bringing from mass eigenfields basis 
to the primed basis and a second transformation U£ that, as already seen, brings from 
the primed basis to the weak basis where the charged lepton mass matrix is diagonal, 
explicitly 

U ZU Ue. (78) 


A generic unitary matrix would be described by nine parameters. However, in the 
specific case of the leptonic mixing matrix, three phases are unphysical since they can be 


absorbed in the charged lepton fields without any observable physical consequence. In this 
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way the leptonic mixing matrix can be parameterised in terms of 6 parameters, 3 mixing 
angles 012,013, 023 and 3 CP violating phases ô, a1, a2. A standard parametrisation is 


then given by 


id 


C12 C13 $12 C13 S13€ 
U = — 819 C23 — C12 S23 S13 €? C12 C23 — $12 S23 513€? S23 C13 - diag (E 1,e'7) , 
S12 523 — C12 C23 $13 ei? — C12 $23 — $12 C23 $13 ei? C23 C13 
(79) 
where s;; = sin6;; and ci; = cos@;. This parametrisation is the same adopted also 


for the CKM matrix, describing mixing in the quark sector, but with the addition of 
the two Majorana phases p and o that cannot be reabsorbed in the redefinition of the 
neutrino fields in the more general case when these are Majorana fields. They cannot have 
any physical effect in low energy neutrino experiments based on kinematical processes 
since these are the same for Dirac and Majorana neutrinos. In particular Majorana 
phases cancel out in neutrino oscillation probabilities and, therefore, cannot be tested 
in neutrino oscillation experiments. At the moment, the only known phenomenology 
that can effectively give us information on Majorana phases is neutrinoless double beta 
(0vG3) decay but since no signal has been found so far, not only we have no experimental 
constraints on the Majorana phases but we do not even know whether they are physical, 
in the case of Majorana neutrinos, or unphysical, in the case of Dirac neutrinos. 
Neutrino oscillation probabilities depend on neutrino mass squared differences and in 
particular they would all vanish if the three neutrinos would be exactly mass degenerate. 
Therefore, the discovery of neutrino oscillations implies that neutrino are massive and 
that at least two of them are non-vanishing and different from each other. Experimental 
results are well described in terms of two mass squared differences corresponding to the 
case of three neutrino masses and do not find evidence of a fourth mass eigenstate. A first 
mass squared difference, Am? = m3 — mî, is measured in solar neutrino experiments 
and in the KamLAND long baseline reactor neutrino experiment. Since matter effects 
in the Sun depend on its sign, this is measured and found to be positive. A second 
mass squared difference was initially measured in atmospheric neutrinos and now also in 
long baseline neutrino experiments. In absolute value it is much larger than the solar 
difference, but the sign is still unknown. For this reason current neutrino oscillation 


experiments do not solve an ambiguity between two possible options: one can have either 


2 


im Z m3 -mÊ > 0 implying mı < mz < ms, or inverted ordering, 


as the 


normal ordering for Am 


for Am? = ms — m% < 0, implying mg < mı < me, where we defined Am?,,, 


mass squared difference with largest absolute value in both cases. 


Al 


The most recent global analyses find for normal ordering |122 


Oi. = 33.44° +0.75° = [31.3°, 35.9°], (80) 
Oo3 = 49.20193 = [40.1°,51.7°], 
O13 = 8.57°+0.12° = [8.20°,8.93°], 
6 = 163° t 04040"), 


sol — 


Amin = (2-517 + 0.027) x 10-9 eV? => Mam = V Amn = (50.2 + 0.2) meV, 


Am? = (7.42 + 0.20) x 10° eV? > my =1/Am2, = (8.6+0.1) meV, 


and for inverted ordering 


te: = BAe SL oes (81) 
O23 = 49.3°99) = [40.3°, 51.8°], 
613 = 8.60°+0.12° = [8.24°, 8.96°], 
ô = —79°726> = [-167°,+8°], 


Am? = (7.42 + 0.20) x 10-3 eV? => Mma = 4/ Am?2, = (8.6 £ 0.1) meV, 


sol — 


Am? = (—2.498 + 0.028) x 10-8 eV? => Mam = V—-AMŽm = (50.00 + 0.25) meV , 


atm 


where we indicated the best fit values with the lo errors and for the three mixing angles 
and ô the 3c ranges as well. 

The measurement of two different neutrino mass scales, the atmospheric and the solar 
ones, implies that we can parametrise the three neutrino masses in terms of just one 
unknown parameter. This can be conveniently chosen as the lightest neutrino mass that 


we denote by my, where 1’ = 1 for normal ordering and 1’ = 3 for inverted ordering. We 


also denote by m3 the heaviest neutrino mass (3’ = 3 for normal ordering and 3’ = 2 
for normal ordering) and with mz the intermediate neutrino mass (2’ = 2 for normal 
ordering and ? = 1 for inverted ordering). If my >> matm then necessarily al three 


neutrino masses become quasi-degenerate. In the case of normal ordering, the hierarchical 


limit is obtained for my = m, < Mg. In this case in the limit mı — 0 one has 
M2 = Mso and M3 = Matm. In the case of inverted ordering, the hierarchical limit 
is obtained for my = m3 < Matm and in this case one has my = mz = Mma, and 
my = mı = Men > m2 ~ Matm [1 — m2,,/ (2 ise) | 


In September 2019 the KATRIN collaboration released first results placing the world 
best direct neutrino mass upper bound on the electron neutrino mass [123]. This has been 


recently improved combining first data taking with second data taking results, obtaining 
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for the first time a sub-eV upper bound |124|, namely 


my, =, [X Veil? m? < 0.8 eV (90%C.L.), (82) 


translating into the same upper bound on the lightest neutrino mass my. If neutrinos 
are Majorana particles then 0vG6 decays are possible. Experiments are sensitive to the 


effective neutrinoless double beta decay neutrino mass defined as 
Mee = |Mizes| = |m UZ, + m2 UZ + m U2, . (83) 


Since a signal was not found so far, they place an upper bound. The most stringent one 
comes from the KamLAND-ZEN experiment: Mee < 165 meV (90% C.L.) [125]. 

Finally, cosmological observations place a stringent upper bound, $`; m; < 0.15 (0.17) eV 
(95% C.L.) [126], on the sum of neutrino masses for normally (inverted) ordered neutrino 


masses, translating into an upper bound 
my <42meV (95% C.L.), (84) 


on the lightest neutrino mass, valid both for normal and inverted ordering. 


4.2 Minimally extended standard model 


We described neutrino mixing assuming that neutrino are massive but without specifying 
how they acquire a mass. For this, one needs a model of neutrino masses. In the SM 
mass terms arise from the Yukawa couplings with the Higgs field 6 = (d+, 4°)? after 
electroweak spontaneous symmetry breaking. We have seen that in the case of charged 
leptons, the Yukawa coupling term Eq. generates the mass term Eq. (64). One can 
then simply think to augment the SM adding RH neutrinos and a Yukawa interaction 
term analogously to the charged leptons. In this way, in a generic lepton flavour basis, 


the total Lagrangian would become £L = L°™” + LY, where 
-Lý = Ly hyp vry ® + h.c., (85) 


with the dual Higgs field defined as ® = icz %*. After electroweak spontaneous symmetry 


breaking, a neutrino Dirac mass matrix is generated by the Higgs vev v 
Mp inj! =U hing (86) 
so that the neutrino Dirac mass term in the Lagrangian can be written as 


A eww 
=L Dismas = Vw Mp yt YRI + h.c.. (87) 
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As for the charged leptons, this can be diagonalised by means of a bi-unitary transformation 
(mathematically, it corresponds to its singular value decomposition) 

muy = (Vit Dan UR) (88) 
where Dmp = diag(mpi, mp2, Mp3). In this way one can switch from the generic primed 
basis of LH neutrino fields vz to the basis of LH mass eigenfields vz; with the unitary 
transformation Vri = Vř;y Vry and, analogously, for the RH neutrino fields by means 
of var = Ug, y Vry. Since the neutrino mass basis coincides with the neutrino Yukawa 
basis where m% is diagonal, the neutrino masses are then simply given by (m1,m2,m3) = 
(mpi, Mp2, Mp3). The leptonic mixing matrix relating the LH mass eigenfields to the LH 
weak eigenfields, considering also the transformation from the primed to the weak basis 
Eq. (74), is then very simply given by U = U£ vy or, in components, Uai = Uf xj! Lir 
If we compare this with Eq. (78), we can make the identification U7 = vi, If charged 
lepton and neutrino Dirac mass terms were diagonal in the same LH neutrino basis then 
Uy =U al and U = I and there would be no neutrino mixing. In the minimally extended 
standard model the mixing is then arising from a mismatch between the charged lepton 
and neutrino Yukawa bases. It is also important to notice that though flavour lepton 
numbers are now not conserved, as observed in neutrino oscillations, the total lepton 
number is still conserved perturbatively. 

The attractive features of this model is clearly that it is minimal and the addition 
of RH neutrinos can somehow appear as a straightforward extension of the SM treating 
neutrinos on the same ground as the other massive fermions. On the other hand, the 
RH components of the Dirac neutrino fields, though they share the same mass with the 
LH components, are sterile degrees of freedom since they are weak interactions singlets. 
Moreover the model presents the following limitations that suggest the necessity of a 


further extension: 


e Though neutrinos also have their RH components, one still requires a rather special 


ad hoc description for neutrino masses, specifically: 


- The lightness of neutrinos compared to all other massive fermions is the result 
of assuming a tiny Yukawa couplings compared to those of the other massive 
fermions. For example, the largest one, explaining the atmospheric neutrino 


mass scale, Matm ~ 0.leV, would be h5 ~ 107:2. 


- The leptonic mixing matrix is the analogue of the CKM matrix in the quark 


sector and, therefore, the observed largeness of neutrino mixing angles, in 
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contrast with the smallness of quark mixing angles, also has to be assumed 


without justification. 


e The model does not provide a solution to the cosmological puzzles and in particular 


to the origin of matter in the universe, in particular: 


- It does not provide a candidate of cold dark matter | 


- There is no mechanism of baryogenesis provided by the minimally extended 
standard model and electroweak baryogenesis would fail as in the SM, since 


the Higgs sector is untouched. 


In addition to these limitations there is also the important observation that, having added 
RH neutrinos, nothing forbids the introduction of a Majorana mass term for the RH 
neutrinos. A possibility would be to impose lepton number conservation, at least at the 
perturbative level but this would require to go beyond the minimally extended standard 
model we discussed, a possibility recently pursued in [127]. The model also contains a 
candidate for dark matter but not a clear solution to the matter-antimatter asymmetry 
puzzle. More generally, one should say that Dirac neutrinos can be still motivated within 
models involving additional ingredients compared to the minimally extended standard 
model circumventing at least some of the limitations we listed but it should be clear 
that in this case the model would not be minimal any more. For example, ways to 
justify the lightness of Dirac neutrinos are found within frameworks with large or 
warped extra-dimensions. Also a mechanism of leptogenesis with Dirac neutrinos 
has been proposed but it still requires some external source, e.g., GUT baryogenesis, 
for the generation of an initial B + L asymmetry so that it anwyay requires an extension 
of the minimally extended sMP4] 

In the following we will consider models where neutrinos are Majorana particles since, 
as it will be clear, they more easily (and successfully) address the problem of the origin 


of matter in the universe. 


?3RH neutrino degrees of freedom can potentially be produced from spin-flip processes and behave as 
unwanted extra radiation at recombination and hot dark matter afterwards. However, with the current 
upper bound on neutrino masses, spin-flip processes are suppressed and the produced abundance of RH 


neutrinos would be completely negligible. 
?4The presence of RH neutrinos can act in a way that part of such initial B—L asymmetry can survive 


sphaleron wash-out. In this way one can have successful baryogenesis models even with B — L = 0. 
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4.3 Minimal (type I) seesaw mechanism 


As we have seen, the mere introduction of RH neutrinos and a Dirac mass term does not 
provide a scenario for the origin of matter in the universe. However, if one abandones 
lepton number conservation at the classical level and even B—L conservation, introducing, 
in addition to the Yukawa coupling term, also a bare (not generated by radiative corrections) 


right-right Majorana mass term M, one has, in addition to the SM Lagrangian, the term 
-== a Te 
= vam = Li hyp vey ® + 5 Yar Mrrvry + he.. (89) 


Now small Latin indexes label LH fields and run from 1 to 3, while capital Latin indexes 
label RH neutrinos and run from 1 to N and, as before, the primed notation indicates 
that we are working in a generic flavour basis. The number of RH neutrinos N should be 
regarded as a free parameter within this general phenomenological framework, since there 
are no model independent theoretical arguments that can determine it. In particular, since 
RH neutrinos are SM singlets and do not carry a gauge anomaly, anomaly cancellation 
requirement does not apply. The Yukawa coupling matrix is then in general a (3, N) 
matrix while the (right-right) Majorana mass matrix is a (N, N) matrix and notice that, 
since RH neutrinos are SM singlets, M can be also a singlet and the term is renormalizable. 

The Yukawa interaction defines the quantum number of the RH neutrinos. Assigning 
L = +1 to all vp fields, the Dirac mass term conserves the total lepton number though not 
the individual lepton numbers or flavours. However, though RH neutrinos carry a global 
charge, they do not carry gauge quantum numbers (colour, weak isospin, hypercharge). 
This is why they can also have a Majorana mass term. However, this now breaks L by 
two units: therefore, the introduction of a Majorana mass term leads to lepton number 
non conservation at the classical level, with processes potentially observable even at low 
energy such as 0vG6 decay. It will also prove important to notice, when we will discuss 
leptogenesis in the next section, that it also clearly violates B — L. 

After spontaneous electroweak symmetry breaking, the Yukawa coupling term produces 
the usual Dirac neutrino mass matrix mh = v h”. However, this time this does not coincide 
in general with the light neutrino mass matrix m,, unless trivially W = 0. Therefore, after 
spontaneous symmetry breaking, we have that the neutrino mass term can be written in 


a generic flavour basis (we omit primed Latin indeces) as 
ieee 
—LY, = TMh VR + 5 RM vet he.. (90) 


Using T; Mp VR = vem vf, we can rewrite the total neutrino mass term as 


pe — 0 mp Ve 
fn = 5 om ( M ) ki 


46 


+ h.c.. (91) 


Introducing the 3 + N RH component neutrino field n’ a = (v, vpr)”, the total neutrino 


mass term can then be written even in a more compact way as 


1— 
-Ly = = ng Mi np t h.c., (92) 


showing that this is in the form of a Majorana mass term also considering that M, is a 


symmetric (3+ N,3 + N) matrix, explicitly 


mu ( al (93) 


Therefore, it can be Takagi diagonalised and the mass eigenstates will be, in general, 
Majorana states] Notice of course that for M = 0 one recovers the Dirac case. The 
most interesting case is the seesaw limit, for M >> mp. One can find a (3 + N,3 + N) 


unitary matrix U diagonalising M so that, performing the transformation 


NR > Np =UNR (94) 
with n% = (v$, Ng)”, one has 
Vy l Ae 
—LY = 5 PR Dm, nR + h.c., (95) 
where Dm, = diag(m1,m2,m3, Mı, ..., Myn). The diagonalisation of M, can proceed in 


two steps. One can first find a block-diagonilising matrix Ug such that 


M, ~ Up DE, Us, (96) 


Di, = ( han (97) 


is a block diagonal matrix and m, is the (3,3) light neutrino mass matrix. The block 


where 


diagonalising matrix Upg is found expanding in powers of £ = mp M7! up to terms O(€?) 


7 1 iE ET —g* 


and at the same time one finds for the light neutrino mass matrix the seesaw formula [131] 


M, = -mp a Mbp. (99) 


25This result is valid in general, not just in the seesaw limit. 
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The lightness of ordinary neutrinos is now explained as an algebraic byproduct deriving 
from the interplay between the electroweak and the new Majorana mass scales. In this 
way the mass spectrum splits into a set of light neutrino masses and into a set of heavy 
neutrinos and correspondingly in the new basis the expression for the Lagrangian 
gets decomposed into a light and a heavy mass term 


| 
Demy vg +- NEM Np. (100) 


1 
=N SZ 
si 2 2 


As a second step, since the light neutrino mass matrix m, is a symmetric matrix, it can 


be further diagonalised by a unitary matrix U% PS]in a way that 
Dm = -U til, (101) 
where Dm = (m1, m2, m3). In this way, using the seesaw formula, we can write 
diag(mı, mz, m3) = U mp a mn UK”. (102) 
The leptonic mixing matrix U is then given by the product 
US. (103) 


If one starts from the basis where the charged lepton mass matrix is diagonal (the weak 


basis) then U£ = I and in this case one simply has U = UY, so that 


1 
diag(mı, m2, M3) = U' mp qm U*. (104) 


As a final step, one can also Takagi diagonalise M with a transformation of the RH 


neutrino fields vgy = U} URI 
Mr y = Upp, Duta URD , (105) 


where Dy = diag(Mr, Mrr, ..., Myn) and we denoted by unprimed capital Latin indexes 
the basis where the Majorana mass is diagonal coinciding with the basis of heavy neutrino 
mass eigenfields. 

The correcting term œ €* ¿T in the diagonal entries in Eq. produces a deviation 
from non-unitarity of the matrix U7, and consequently of the leptonic mixing matrix U. 
Therefore, within the type-I seesaw mechanism, one can potentially have non-unitarity 


effects in neutrino oscillations within a standard three neutrino mixing approach. However, 


26 As discussed in detail in footnote 3, mathematically this is called Takagi diagonalisation of a complex 


symmetric matrix, different from the usual eigenvalue decomposition. 
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in high-scale seesaw models these effects are highly suppressed but they can show up in 
seesaw models at the TeV scale |132H134|. 
It is then finally useful to write the light and heavy mass eigenfields, that are Majorana 


fields, as a linear combination of the LH and RH neutrinos: 


Vi 


Vin + Ver = wee + (Ulva)! = [Viy Ei y VRI + (Vii ory vrr )'] 5 (106) 
Nr = Nirt+ Nip = Uae vrr + Uri ver) | + [Uri Shi gt vry + Ur Enp Vij)" | : 


This clearly shows how the light neutrinos are dominantly LH neutrinos but also have a 
subdominant RH neutrino component and vice-versa heavy neutrinos are dominantly RH 
neutrinos but also have a small subdominant LH neutrino component and the mixing is 
described by €. For high-scale seesaw models € is too tiny and the mixing does not give 
any sizeable phenomenological effects. However, for low-scale seesaw models the mixing 
becomes larger and this has a host of potentially testable physical applications ?"| 

First of all the heavy neutrinos are not any more fully sterile and they can couple 
to SM particles and can be potentially detected in colliders. Since so far no signal was 
found, upper bounds on the mixing parameters that describe this active-sterile neutrino 
mixing are placed. However, cosmologically, even small mixing angles can be effective 
in producing a heavy neutrino that can be long-lived and play the role of dark matter 


particle. We will discuss this cosmological application in Section [6] 


4.4 Two extreme limits of type-I seesaw models 


The seesaw mechanism by itself does not specify the masses of the RH neutrinos. Model 
independently, these have to be regarded as free parameters and, in general, cannot be 
expressed in terms of the low energy neutrino parameters. This interesting possibility 
would then require some assumption on the neutrino Dirac mass matrix. If one considers, 
as an illustrative example, a one generation toy model, the seesaw formula reduces simply 
to M = m?,/m,. If one assumes mp ~ Mpgw ~ 100 GeV, then one obtains M ~ 1015 GeV 
using My ~ Mso ~ 10meV. This can be regarded as an encouraging exercise since the 
light neutrino mass scales would be obtained from the electroweak scale and a scale very 
close to the grand-unified scale. However, it is clearly not compelling and lowering the 
neutrino Dirac mass scale below the electroweak scale one can lower the RH neutrino 


mass scale by many orders of magnitudeļ”] 


2T For a general discussion on the phenomenology of RH neutrinos see i 
28 As we will see, cosmologically, it will be particularly interesting to consider mp ~ 1eV, my ~ 1074eV 


implying M ~ 10keV, though clearly in this case the neutrino Dirac neutrino mass, and therefore the 


related Yukawa coupling, is again much smaller than in the case of all other massive fermions. 
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Accounting for a realistic three generation case, the indetermination on the RH neutrino 
mass spectrum is even more evident since the neutrino Dirac mass matrix contains fifteen 
parameters while, on the other hand, there are only nine low energy neutrino parameters 
in m, and we do not even test all of them considering that 0v6 decay experiments give 
information on just one quantity involving two Majorana phases. 

However, using the seesaw limit and the seesaw formula Eq. (104), there are some 
interesting limits, corresponding to definite theoretical assumptions, where it is easy to 
extract simple relations expressing the heavy neutrino masses in terms of the low energy 


neutrino parameters. 


4.4.1 Degenerate Dirac neutrino mass spectrum 


Let us for definiteness consider the most attractive case N = 3, motivated in different 
models. A very simple assumption is to assume mp = AJ, where A is an overall mass scale 
and J is the identity matrix in flavour space. This form is highly symmetric and invariant 
in any basis and emerges quite straightforwardly from symmetry flavour models (136). 
In this case one can immediately invert the diagonalised see-saw formula Eq. (104), 
since the leptonic mixing matrix U cancels out, finding for the RH neutrino masses 
eee, Be ik= (107) 
m3 m2 Mı 
simply the inverse of the light neutrino masses. In Figure 6 we show a plot of the three 
heavy neutrino masses versus the unknown lightest neutrino mass m; in the case of normal 


ordering. 


4.4.2 SO(10)-inspired models 


Another important example where it is possible to express the RH neutrino mass spectrum 
in terms of the low energy neutrino parameters is the case of SO(10) inspired models [137]. 
In this case, in contrast with the previous limit case, one assumes that the neutrino Dirac 
mass spectrum is very hierarchical and more precisely the neutrino Dirac masses are of 


the same order of magnitude as the up quark masses, explicitly 
Mp1 = Q1 Mup, Mp2 = Q2 Mcharm, MDZ = Q3 Miop ; (108) 


where the three a; are O(1) factors. This implies approximately mp; : mp, : mp3 ~ 
Wigs : Mecharm : Mtop ~ 1 : 1025 : 10°. With this assumption, one can derive the following 


expressions for the three RH neutrino masses [138}/139| 


[Myl =1 


M; ~ mhs |(™,") 33] , (109) 


2 
M, ~ Mp1 M ~ MDə < , 
Mı Mma M3 |(Mz+)33] 
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Figure 6: Spectrum of RH neutrino masses versus m; in a discrete flavour symmetry model 
where mp = AJ in the case of normal ordering (from |136]). In the plot the parameter 
ys = A/v. The vertical band indicates the allowed range of values for m; in the specific 
A4 model considered in [136]. This fixes quite precisely the RH neutrino masses. 


where m, is the light neutrino mass in the neutrino Yukawa basis, where mp is diagonal. 
Notice that these expressions are not valid in the vicinity of m,1, = 0 or (m;')33 = 0 
or both. In these cases the spectrum two or even all three RH neutrino masses become 
quasi-degenerate. The expressions we gave work well as far as the RH neutrino masses 
are hierarchical. In the singular points there are so called crossing level solutions [140]. 
However, it has been noticed that in these cases the light neutrino masses become highly 
fine-tuned and, therefore, they should be considered as some very special cases. We will 
be back on this point when we discuss leptogenesis. 

If we consider the neutrino Dirac mass matrix in the weak basis, where the charged 
lepton mass matrix is diagonal, and diagonalise it by means of the bi-unitary transformation 
in Eq. (88), we can write (a = e, u, T; I=1,2,3; k = 1, 2,3) 


MDal = Vis Dmpk URkJ « (110) 


The unitary matrix V, now transforms the LH neutrino fields from the weak to the 
neutrino Yukawa basis, where the neutrino Dirac mass matrix is diagonal. It also transforms 


the light neutrino mass matrix from the weak to the Yukawa basis, explicitly 


My = Vrm, VF. (111) 
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1074 1078 107? 107 1 
my, (eV) 
Figure 7: Typical RH neutrino mass spectrum in SO(10)-inspired models for V; = J, 


Qy Q2 Q3 1, normal ordering, best fit values for the three mixing angles (see 
Eq. (80)) and p = ø = ô = 0. The solid lines are the numerical solutions and the dashed 


lines are the analytical expression given in the text. 


The matrix Vz; has to be regarded as the analogous in the leptonic sector of the CKM 
matrix Vox in the quark sector, bringing from the basis where down quark mass matrix 
is diagonal to the basis where the up quark mass matrix is diagonal. In SO(10)-inspired 
models it is then assumed that the matrix Ve ~ Voxyy. Since the largest mixing angle 
in Vcxm is the Cabibbo angle 0e ~ 13°, the dependence of the RH neutrino masses 
on Vz is quite mild and with a good approximation one can simply use Vsp ~ J. In 


Fig. |7| we show a typical example of RH neutrino mass spectrum in SO(10)-inspired 


models for V; = J and a, a2 a3 1. In particular, the three M7; are plotted 
as a function of the lightest neutrino mass mı in the case of normal ordering. One can 
notice the very strong hierarchical RH neutrino mass spectrum that is obtained. It is also 
useful, from the analytical expressions Eq. within the approximation Vz; = J and 
for a, = a2 = a3 = 1, to derive the RH neutrino mass spectrum in the hierarchical limit 


for mı — 0 


2 2 
m m m 
Mı = DL , Mz = — 2? >, M; ~ 2 533 Sia, (112) 
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and in the quasi-degenerate limit for m1 >> Matm (implying Mı ~ M ~ m3) 


2 
mMm p — 
M ~ i ta 1], (113) 
Mı 
2 
m : 
Mə ~ D2 1 2 2ip 1 
2 ma | + cigle )|, 
2 
m ; 
M; ~ = |i + s3 (e — 1) + l. 


Mı 


It can be noticed how in the hierarchical limit one has M3 œ mj", a necessary consequence 
of the seesaw formula. In this limit the heaviest RH neutrino decouples from the seesaw 
formula and effectively one obtains the two RH neutrino case with a considerable reduction 
in the number of seesaw parameters, precisely from eighteen to eleven (141). For example, 
for mı = 10-*eV one can see from the plot that Mı ~ 10°GeV, Mə ~ 10! GeV and 
M ~ 10'© GeV. We will see in the next section that such a RH neutrino mass spectrum 
will have important implication for the calculation of the matter-antimatter asymmetry 
within leptogenesis. Here we just notice that M3 ~ Meur and this is indeed the most 
straightforward expectation for the RH neutrino mass spectrum from a grand-unified 


model such as SO(10). 


4.5 Beyond the minimal seesaw mechanism 


Some of the results we derived within the type-I seesaw mechanism can be generalised 
within an effective theory approach pioneered by Weinberg [142]. In particular, the 
lightness of the neutrino masses from the seesaw formula can be understood in a more 
general way. 

We have seen that starting fro the type-I seesaw Lagrangian (90), after the diagonalisation 
of the total mass term, we could express the mass terms in the mass eigenstate basis 
finding that the light neutrino mass term is given by a left-left Majorana mass term (see 


Eq. (100)) 


-Lp = Vimy VF. (114) 


This result is more general. First of all one could think to introduce directly a left-left 
Majorana mass term of this kind without introducing RH neutrinos and a Dirac mass 
term. This would seem a very economical way not requiring apparently any extension of 
the SM. However, things are not that simple. First, in that case m, would have still to be 
fine-tuned to reproduce the tiny neutrino masses. Even more importantly, the Lagrangian 
has to be not only Lorentz invariant but, within the SM, also SU(2); x U(1)y invariant 
not to spoil the SM renormalisability and the only way for the mass term in Eq. to 
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be gauge invariant is to have m, transforming as a SU(2),, triplet with Y = a] However, 
the SM does not contain any weak isospin triplet with Y = 2. For this reason it is not 
possible to build, within the SM, a renormalizable Lagrangian term which can generate 
a left-left Majorana neutrino mass. Therefore, in order to have a gauge invariant term 
without spoiling renormalizability, one necessarily needs to extend the SM introducing 
new fields. We have seen the solution provided by the type-I seesaw mechanism. Another 
possibility is to introduce an isotriplet Higgs field that couples to the lepton doublets 
through Yukawa couplings ha. After spontaneous symmetry breaking such a Yukawa 
term would indeed generate a left-left Majorana mass term. 

However, there is a more general way, based on effective theory, to generate m, and also 
elegantly predicting the lightness of neutrino masses from a generalised seesaw formula. 
The left-left Majorana mass term can be indeed regarded as generated by the five 
dimensional (lepton number violating) Weinberg operator |142], 


f= (Le &*) (Gt Lg) + h.c., (115) 


where gag is a symmetric dimensionless coupling constant matrix, and A is the energy 
scale of new physics (e.g., the mass of a new super heavy particle). This is the lowest 
dimensional operator one can construct with the SM fields respecting SM symmetries. 
After spontaneous symmetry breaking this operator does indeed generate a left-left neutrino 


Majorana mass 
m, — I vV mb 
vaß — a A 
A A 


This is quite interesting since if A > mp one can naturally suppress the neutrino masses 


(116) 


at the level we measure, and in this way the lightness of neutrino mass is associated to 
the existence of new physics at very high energy scale: Eq. can then be regarded 
as a generalised seesaw formula. Since the Weinberg operator is non renormalizable, it 
has to be regarded as an effective operator generated at energies much below the scale 


of new physics A. There are only three ways to generate the Weinberg operator at tree 


level [144]: 


1. Type-I seesaw mechanism: with a SU (2) singlet fermion as a mediator [131]; 
2. Type-II seesaw mechanism: with a triplet scalar as a mediator [145]; 


3. Type-III seesaw mechanism: with a triplet fermion as a mediator [146]. 


29A neutrino LH chiral field vz has isospin third component T3 = 1/2 and hypercharge Y = —1 (see 
for example Table 3.5 in 143]). This implies that vg vf has Jz = 1 and Y = —2. 
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If one assumes the couplings gag ~ 1, then necessarily A ~ Mgur in order to reproduce 
the solar and atmospheric neutrino mass scales. This is obviously quite an attractive case. 
The drawback is that such high scale of new physics, and corresponding high-scale seesaw 
models, escapes direct detection at colliders. In the last years there was intense activity 
in considering alternative possibilities where neutrino masses might be generated by new 
physics at the A ~ TeV scale or even below, directly accessible at colliders such as the 


LHC. In this case there are three main strategies: 


(i) Low-scale seesaw models: couplings gag are very small; their suppression is then 
usually justified in terms of small lepton number violation (e.g., linear and inverse 
seesaw mechanisms, see for example [132}) or it is not justified, small couplings are 


just imposed by hand; 


(ii) Neutrino masses are generated radiatively and their suppression is guaranteed by a 


combination of loop integrals and electroweak scale masses entering the diagrams 


(see for a recent review). 


(iii) Neutrino masses from the Weinberg operator are forbidden (gag = 0) but they are 
generated by higher dimensional operators|>)| 


We have then a picture with a very broad variety of possibilities. We will now discuss how 
the necessity to explain the origin of matter in the universe can shed light on which of 
these options has more attractive features and can be regarded as more successfully. Here 
we just want to comment that low scale seesaw models are certainly very attractive from 
a phenomenological point of view but on the other hand the null results from searches 
of new physics at colliders and with many other different experimental strategies so far 
make these models more and more constrained experimentally. As we will see, cosmology 


will draw us to similar conclusions. 


5 Leptogenesis 


Neutrino physics offers an elegant solution to the problem of the matter-antimatter 
asymmetry of the universe. Moreover, among the long list of baryogenesis models, has 
the clear advantage that relies on the only new physics that has been experimentally 
established: neutrino masses and mixing. Since the original proposal [29], many variants 


have been proposed. Here we mainly focus on the minimal model of leptogenesis since 


3°This case has common features and overlap with (ii) |147]. 
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it remains the most attractive scenario, supported, as we will point out, by a few simple 
experimental facts and also because non-minimal extensions, usually low-energy scale 
scenarios, are not supported by the lack of new physics in experiments (beyond neutrino 
oscillations) so far. Moreover in the minimal scenario the measured value of the baryon-to-photon 
number ratio points to a RH neutrino mass scale, M ~ 101071 GeV, that emerges quite 


naturally in models. 


5.1 The minimal scenario of leptogenesis 


The minimal scenario of leptogenesis PJ is based on two assumptions, one on the particle 


physics side and one on the cosmological side: 
i) Type-I seesaw extension of the SM; 
ii) Thermal production of the heavy neutrinos (thermal leptogenesis). 


As discussed, the explanation of the lightness of neutrinos is well addressed by the 
type-I seesaw mechanism that predicts the existence of heavy, dominantly RH, neutrinos. 
High-scale seesaw models are difficult to test since the parameters associated to the heavy 
neutrinos, like their masses, escape the experimental investigation. However, the heavy 
neutrinos, though typically decaying prior to the big bang nucleosynthesis onset without 
affecting standard cosmology, can leave a crucial imprint at present: the matter-antimatter 
asymmetry of the universe in the form of the observed baryon asymmetry. 

Leptogenesis is another important example, like GUT baryogenesis, of baryogenesis 
from heavy particle decays[?] Their decay can indeed violate CP since in general the rate of 
decays into leptons T(N; —> Lr+¢ġ') can be different from the decay rate into anti-leptons, 
T(N > Lr +¢), in way that each N; decay will produce, on average, a B — L asymmetry 
in the form of a lepton asymmetry given by the the total CP asymmetry, defined as 


[(N; > Lr + ¢') -T(N; > Lr 4+ 4) 
T(N; +> L,+ 61) +T(N; + Lr+ o) l 


E= — (117) 
Compared with the general definition Eq. (i), we now have sign(Ag_;) = —1. As we 
will see, when flavour effects are included there is an additional source of CP violation. 
For the time being we neglect flavour effects, we will see how these can change the results. 

Even though the B — L asymmetry is initially injected in the form of a lepton 
asymmetry, if the reheat temperature of the universe, the initial temperature at the onset 


31For dedicated reviews see f 
et 


32We can then here specialise the general formalism and results discussed in Section 2.4. 
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of the radiation dominated regime after the end of inflation, is above the out-of-equilibrium 
temperature of sphalerons, this will be redistributed among all SM species and in particular 
also among quarks, so that a baryon asymmetry is also generated. Chemical equilibrium 
150| ultimately enforces L œ (aspn— 1) (B—L) and B ~ aspa (B—L), with aspn = 28/79 ~ 
1/3. 

Notice that also leptogenesis, likewise electroweak baryogenesis, crucially relies on 
sphalerons to generate a baryon asymmetry. In both cases the initial injected asymmetry 
is not a baryon asymmetry but a chiral asymmetry in the case of electroweak baryogenesis 
and a lepton asymmetry in the case of leptogenesis|**] The important difference is that in 
the case of leptogenesis B — L is violated by the decays of the RH neutrinos, while it is 
conserved in electroweak baryogenesis. In this respect leptogenesis is directly relying on 
new physics while in the case of electroweak baryogenesis this could potentially be viable 
even in the SM. However, as we discussed, electroweak baryogenesis fails to a quantitative 
level because the experimental results on the Higgs mass and CP violation in quarks show 
that departure from thermal equilibrium and CP violation are not enough and for this 
reason new physics is needed. 

The final B — L abundance in a portion of comoving volume will be then given by the 


sum of the contributions from all RH neutrino species 


1N 
Net =) erri, (118) 
T 


where the «{’s are the final efficiency factors taking into account the wash-out from 
inverse processes, the efficiency in the RH neutrino thermal production and also contain 
geometrical factors depending on the flavour structure. For definiteness we will consider 
the case N = 3, also because this is the most attractive case emerging in different models. 


From the final B — L asymmetry one can then calculate the baryon-to-photon ratio 


at present 
Ni 
lep B-L 
= Asph —— - 119 
1B0 ph N-o (119) 


If we normalise the abundances in a portion of comoving volume containing one RH 
neutrino in ultrarelativistic thermal equilibrium, prior to the onset of leptogenesis, at 
some fiducial initial time t; such that Ti > Mr, where T; = T(ti) and My, is the mass 
of the heaviest RH neutrino contributing to the final asymmetry, then the number of 


photons at this stage is simply Ni = 4/3. Assuming a SM particle content after RH 


33For some reason this aspect is always emphasized for leptogenesis but not for electroweak baryogenesis 


that after all could be also called ‘chiralogenesis’. 
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neutrino decays and a standard cosmological expansion, implying entropy conservation, 
then one has No = Ni g4 /gso, Where the number of (entropy) degrees of freedom at 
T > M; is given by the number of SM degrees of freedom plus the contribution from the 
RH neutrino degrees of freedom, simply gb = gsm = 427/4 +7/4 = 108.5 and the number 
of entropy degrees of freedom at the present time is gso = 43/11 (151). In this way, one 
finds Nyo = 4774/129 ~ 37. Finally, with the normalisation of abundances discussed and 
adopted in Section 2.4, one finds 
ep 28 129 
"Bo = zo * A774 
This is the quantity predicted by leptogenesis that has to match the measured value in 
Eq. (5). Clearly one needs to calculate Nj, and from Eq. one can see that for 
this one needs to calculate the total CP asymmetries €r of the RH neutrinos and their 


Nig = 0.96 x 107° N z. (120) 


efficiency factors. The first depend just on the particle physics and in particular on the 
neutrino Dirac mass matrix and, very interestingly, on the RH neutrino mass spectrum. 
The efficiency factors are a mix of particle physics and cosmology, they clearly require a 
solution of kinetic equations. 

A perturbative calculation from the interference of tree level with one loop self-energy 


and vertex diagrams gives for the total CP asymmetries 155 


3 5 Im [(h' h)fs] E(z3/21) 


ga | (121) 
ce ar Co) Vaser 
where we introduced 
2 1l+2 2-2 
T 1 l — 122 
g(a) = 52 [a +a) m (HE) -E (122) 


and defined zy = M?/M?. An important property of total CP asymmetries is that, 
depending on the combination h' h, the dependence on the leptonic mixing matrix cancels 
out. We will see that this will have an important consequence. 


The total decay rates are given by 
I? =IT(N; > Lr + ¢') +T(N; = £7 4+ o) (123) 
and at zero temperature they correspond to decay widths given by 
rp =IP(f =0) = an (hi h)rr. (124) 


The total decay rate can then be obtained multiplying the decay width by the averaged 


dilation factor, 


rp =T? (=) l (125) 


YI 
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The averaged dilation factor is given by Eq. (21), where now z has to be replaced by 
zy = M,/T. The decay factor, introduced in the general case in 2.4, can now be written 


as Dr =? /(H zr). The decay term can be conveniently expressed in terms of the decay 
parameter (see Eq. (16)) as 


ine. 
Dy(z1) = Kızı (=) = Krz?. (126) 


YI 


The RH neutrino decay parameters are interestingly related to the light neutrino 


masses and can be recast as 


My A Be v (h h)rr 
Kr = — th = —___ 127 
r=, with s (127) 
where the Mpr’s are the effective neutrino masses and 
16 5/2 SM 2 
m, = Z V9 © ~107mev (128) 


3V5 Mp 


is the equilibrium neutrino mass [157]. The seesaw formula, combined with neutrino 
oscillation experiment results, typically favours My in the range Msoi Matm ~ (10-50) meV, 
corresponding to K; ~ 10-50. More specifically, one cannot have more than one decay 
parameter Ky; < Mso/M„ ~ 10. This point plays an important role both in leptogenesis 
and also when we will discuss the possibility that one heavy seesaw neutrino could play 
the role of dark matter. 

As discussed in Section 4, the type-I seesaw extension of the SM introduces eleven 
new parameters in the case of two RH neutrinos and eighteen new parameters in the case 
of three RH neutrinos. Even in the case that Ov decay experiments find a positive 
signal and measure Mee and cosmological observations measure the lightest neutrino 
mass My, this will determine only eight parameters in the low energy neutrino mass 
matrix. The RH neutrino masses and the complex mixing angles in the orthogonal matrix 
will still escape the experimental information from low energy neutrino experiments. 
Leptogenesis can then be regarded as a special very high energy neutrino experiment that 
took place once in the early universe and that is sensitive also to high energy neutrino 
parameters. However, in the end it just introduces one additional experimental constraint, 
the successful leptogenesis condition nso = npo that, in general, is still not sufficient to 
over constraint the seesaw parameter space making model independent predictions, even 
in the case of two RH neutrinos (N = 2). 

To this extent it is necessary to reduce the number of parameters, increasing the 


predictive power, and/or to look for additional phenomenological tests. In the first 
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case one can either introduce some theoretical condition within a model or a class of 
models or make some approximations and reasonable assumptions in the calculation of 
the asymmetry that can lead to some experimental prediction on low energy neutrino 
parameters. In the second case one can either lower the scale of leptogenesis, with or 
without going beyond a type-I seesaw mechanism, and/or try to reproduce some other 
phenomenological observation. The most remarkable ohenomenological observation one 
can try to combine with the matter-antimatter asymmetry is certainly the dark matter of 
the universe. In this way one would obtain a unified model of neutrino masses, leptogenesis 
and dark matter, a model that would solve the problem of the origin of matter in the 
universe. In the following, we will discuss all these different strategies but of course the 
possibility to explain also dark matter will be our main focus. We start from the simplest 
and most traditional approach that leads to a prediction on the absolute neutrino mass 
scale from leptogenesis, based on a set of assumptions and simplifications in the calculation 


of the asymmetry: so-called vanilla leptogenesis [148]. 


5.2 Vanilla leptogenesis 


Vanilla leptogenesis is the simplest scenario of leptogenesis and it is an important example 
of how, just imposing the successful leptogenesis condition, one can obtain predictions on 
physical observables, specifically an upper bound on the absolute neutrino mass scale. It 


relies on the following set of assumptions: 


(i) The flavour composition of leptons and anti-leptons produced in the decays and 


inverse decays of RH neutrinos does not influence the value of the final asymmetry; 
(ii) The RH neutrino mass spectrum is hierarchical and more precisely Mz Z 2 M1; 
(iii) The lightest RH neutrino decay parameter Kı > 1; 
(iv) There are no fine-tuned cancellation in the seesaw formula. 


The first three assumptions imply that the final asymmetry is dominated by the lightest 
RH neutrino (Nj;-leptogenesis), so that the expression (120) gets specialised into 


nish ~ 0.96 x 10772, K8"(Ky, my), (129) 


where notice that the final efficiency factor for the lightest RH neutrino decays depends 
just on the lightest RH neutrino decay parameter and, quite interestingly, on the absolute 


neutrino mass scale expressed in terms of the lightest neutrino mass mı (either mı for 
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normal ordering or m3 for inverted ordering). The lightest RH neutrino efficiency factor 
can be calculated solving a simple set of just two rate equations, one for the N; abundance 


and one for the B — L abundance 
dNy, 


= =), (Ny, = NS 130 
dz it Ny No ( ) 
dNp_ o 
= = = —e Dı (Nm, — NA) — WP Np-i, (131) 
1 


where, analogously to the decay term D1, we defined the inverse decay wash-out term as 


1 
mZ = K z . 132 
Wi 2 H(z) 21 4 1721 Kiı(zı) ( ) 


In the last expression we used that the inverse decay rate is related to the decay rate (see 


Eq. (125)) simply by 


Ny, (z1) 1 
rP (2) SPP (zi) Na) = T? 2 Ko(z)- (133) 
L 


A simple solution for the final efficiency factor can be written in the form 


! W Mı m 
Ki(Kı mi) = «i(K1) exp -2 1010 GeV (Z) ’ (134) 
where 5 
f _ -1z pam] 
kens = (1- B , 135 
Ky ( 1) zp(Ka) Ky e ? ( ) 
zp(Ky) 244 K2% eR (136) 
and 


m = yo. (137) 


The exponential is due to wash-out from AL = 2 scatterings and the constant w is given 
by 
_ 9/5 Mp 1078 GeV? 
4r9/2 gi VI vt 
This value implies that for Mı < 10! GeV the exponent is negligible but for Mı = 


M 


~ 0.186. (138) 


1013 GeV it becomes large, inducing an exponential suppression. Finally, using the assumption 


(iv), one can derive an upper bound on the total CP asymmetry |159 


= 3 Mı Matm = Mı Matm 
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Figure 8: Allowed region in the plane Mı versus mı. The lower bound on the lightest RH 


neutrino mass implies a lower bound on the reheat temperature. 


In this way one has an upper bound on the asymmetry depending just on three parameters: 
Kı, my, M, explicitly: 


nse < 0.96 x 10722 (M1) K® (K1, my). (140) 


Imposing the successful leptogenesis condition, from the upper bound one finds 
the allowed region shown in Fig. |8| in the plane Mı versus mı Pl One can first of all 
notice the lower bound on the lightest RH neutrino mass M; = 10° GeV [68||159], coming 
from the upper bound on £, and also, interestingly, the existence of an upper bound on 
the neutrino masses mı < 0.1eV [64], now confirmed by the cosmological upper bound 
obtained within the ACDM model. This nice confirmation is certainly an encouraging 
result supporting vanilla leptogenesis but it also provide a clear simple example of how, 
under certain assumptions, one can get predictions on experimental quantities. 

Another encouraging feature of vanilla leptogenesis is that the assumption of strong 
washout, Kı > 1, is supported by neutrino oscillation experiments. As we said, the 


seesaw formula typically favours Mı in the range Mgo1-Matm ~ (10-50) meV corresponding 


34The figure is obtained for normal ordering but a very similar result is also found for inverted ordering. 
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to Ky ~ 10-50 #5] For this reason having Kı < 1 is a very special situation | This range 
of values of K, corresponds to efficiency factors kf ~ 1073-107? that still allows successful 
leptogenesis if Mı is above the lower bound. On the other hand, for such values of Kı the 
efficiency factor is independent of the initial N; abundance and moreover any pre-existing 


asymmetry is very efficiently washed-out, since its final value is exponentially suppressed 


Nei = NBA ee, (142) 


Therefore, neutrino oscillation experimental results and the upper bound on neutrino 
masses from cosmology greatly support the vanilla leptogenesis scenario. However, there 
is an unpleasant point in these results: the lower bound on M; rules out SO(10)-inspired 
models we discussed in the previous section where, as we have seen, one has M, ~ 10° GeV 
unless one considers special crossing level solutions [140]. However, these solutions imply 
a large fine-tuning in neutrino masses from the seesaw formula [138}[139]. We will see 
that the inclusion of flavour effects dramatically changes this conclusion opening up 
a new interesting scenario of leptogenesis where the asymmetry is produced from the 
next-to-lightest RH neutrinos (V2-leptogenesis). 

The existence of a lower bound on the lightest RH neutrino mass also translates into a 
lower bound on the reheat temperature of the universe that is slightly more relaxed since 
the asymmetry is produced quite sharply around a temperature Tg ~ Mı/zg ~ M,/(2-5). 
The existence of such a lower bound on the reheat temperature can be a problem in natural 
gravity mediated supersymmetric models because of the gravitino problem whose 
solution typically requires reheat temperatures below 10?GeV. However, the current 
experimental status simply does not support natural supersymmetric models and points 
to a situation where in any case some level of fine-tuning in the stabilization of the 
electroweak scale is unavoidable and the upper bound from the gravitino problem gets 
either relaxed or simply does not hold, in any case it is certainly less compelling than 
initially thought. 

The existence of the lower bounds on M, and on the reheat temperature have motivated 
an intense research on how to evade them departing from the simple vanilla leptogenesis 


scenario. One of the motivation, as we mentioned was the gravitino problem, but in 


3°Tn this range a simple and very good analytical fit for the final efficiency factor is given by: 


(141) 


36Recently it was shown that for randomly uniformly (in flavour space) generated seesaw models 
respecting neutrino oscillation experimental results, the probability to have Kı < 1 is less than 0.1% ; 
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recent years the possibility to have successful leptogenesis at the TeV scale gained a lot 
of attention for two reasons. A first reason is that the presence of RH neutrinos tends to 
destabilise the electroweak scale unless some degree of fine-tuning is requested. Imposing 
that the fine-tuning is not higher than percent level, one obtains the upper bound 
M; < 108 GeV. 

At the same time, the possibility to have leptogenesis at the TeV scale opens the 
opportunity to test the existence of RH neutrinos at colliders. Of course this second 
reason is justified from the first one, essentially naturalness. Since the LHC is placing 
severe lower bounds on the scale of new physics, this has, however, risen doubts on 
naturalness driven arguments and in the last years these two motivations became less 
compelling. 

A first way to lower the scale of leptogenesis is resonant leptogenesis 164]. In this 
case it just relies on relaxing the assumption of a hierarchical RH neutrino spectrum since 
the CP asymmetries get enhanced when the lightest and next-to-lightest RH neutrino, 
or all of them, get quasi-degenerate. However, barring fine-tuning in the type-I seesaw, 
Yukawa couplings get necessarily too small at the TeV scale to make the RH neutrinos 
detectable. Therefore, resonant leptogenesis needs to work in combination with some 
additional extension of the seesaw Lagrangian. One possibility is that RH neutrinos have 
extra gauge interactions and a typical example comes from gauging a Ug_;(1) symmetry 
getting broken at low scale [165]/166]. In this case one gets a rich phenomenology at 
colliders with Z’ detection. Currently stringent lower bounds have been placed at the 
LHC [135]/167}{169]. Another popular way to lower the scale of leptogenesis and at the 
same time to have chances to detect RH neutrinos at the LHC relies on going beyond 
a pure type-I seesaw mechanism, for example considering leptogenesis with both type-3 
and type-1 contributions to the seesaw formula or in recent years leptogenesis within 
linear and inverse seesaw mechanism has been also explored [171]. In this case the one 
can justify large Yukawa couplings that, as we said, in a type-I seesaw mechanism would 
be highly fine-tuned. Finally, another popular low-scale scenario of leptogenesis relies on 
a production of the asymmetry from RH neutrino mixing rather than from decays, the 


so-called ARS scenario |172}174). 
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5.3 Flavour effects 


Within the minimal scenario defined in Section 5.1, the most important generalisation of 
the vanilla leptogenesis scenario comes from considering flavour effects The 
flavour composition of lepton in final states is neglected in vanilla leptogenesis. However, 
in certain situations this can dramatically change the calculation of the final asymmetry. 
If the mass of the decaying RH neutrino is much above 10!? GeV, then active processes 
in the early universe are flavour blind. In this case the unflavoured approximation and 
the simple Boltzmann equation Eq. for the calculation of the asymmetry provides a 
fair description. However, for lighter RH neutrinos with masses in the range 10° GeV < 
Mı « 10!? GeV, Yukawa interactions of the lepton quantum states |¢,;) produced in the 
N,-decays with the RH tauons in the thermal bath, will measure, on average, the tauon 
component before they inverse decay. In this way the inverse decay wash-out will now act 
independently on the tauon component and the electron plus muon coherent component 
of the |@,)’s. One has then to track separately the tauon asymmetry A, = B/3 — L, and 
the electron plus muon asymmetry Ae}, = 2B/3 — Le — L,,. In this two fully flavoured 
regime the set of Boltzmann equations and is then replaced by 


dN. 

A = -Di (Nn, — NB), (143) 
Z 

dNa 


t = Ele Dy (Ny, — Ny) — Na. pie WP , 


c= Eiu D1 (Ny, — Nxt) — Na, Plegu WP - 


Here pł, and Pin are the probabilities, at tree level, that a quantum lepton state |41) 
is measured by tauon Yukawa interactions respectively as a tauon or an electron plus 
muon lepton state and in terms of the neutrino Dirac mass matrix they can be simply 
calculated as pl, = |mpia|?/(m', mp). From these one can also introduce the flavoured 
decay parameters Kia = pł, Kı. In Eqs. the flavoured CP asymmetries £1a are 
defined as _ 
liamli 
G+ 
The flavoured decay rates are related to the total decay rates by lia = pial) and Tia = 


(144) 


Ela 


Pia l1, where pia and pja are the probabilities that the lepton and anti-lepton quantum 
state |¢,) and |¢;) are measured in the flavour a by the interaction with a RH charged 


lepton a. Notice that in general Apia = Pia — Pi, Æ 0 and this introduces an additional 


3TFor a recent review on flavour effects see . Before the the general discussion in , flavour 
effects were specifically considered within resonant leptogenesis and N2-leptogenesis 180]. 
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source of CP violation. For this reason one can have the generation of an asymmetry even 
if the total CP asymmetry vanishes 
Solving the set of Boltzmann equations one finds the final flavoured asymmetries and 


from these the final B — L asymmetry is simply given by their sum 


Nb- = Na. M NA (145) 


ety ` 


Notice that in the limit of vanishing wash-out, for Kı < 1, the sum of the two equations 
for the flavoured asymmetries simply reduces to the Eq. and, therefore, one obtains 
exactly the same solution as in the unflavoured regime. Therefore, flavour effects can have 
an impact on leptogenesis predictions only if wash-out is non negligible, for Kı = 1 [158]. 
In this case one can arrive to the following approximate expression for the final B — L 


asymmetry |181|: 


Apir 
NË, = 26, k (Ky) + A 


[E (Kir) — si (Kieay)] - (146) 
The first term is just twice what one would obtain in the unflavoured case, an effect 
that was already spotted in early works studying flavour effects [67], but the second 
term is something genuinely arising from the account of flavour effects and can generate 
dramatic changes in the result compared to the unflavoured case. The quantities Ap;, and 
APie+p = —Apir represent the difference in the flavour composition of lepton quantum 
states and anti-lepton quantum states CP conjugated and rise from radiative corrections. 
If the flavoured decay parameters in the two flavours the difference in the second term 
is non-vanishing and the Apia can, under certain conditions, be comparable or even 
higher than the total CP asymmetry. In this situation the final B — L asymmetry can 
be strongly enhanced compared to the unflavoured case. In this respect it can be shown 
how CP violating phases, and in particular the Majorana phases, play an important 
role in making in a way that difference between the two final efficiency factors in the 
two flavours does not vanish. Therefore, flavour effects introduce a dependence of the 
final asymmetry on the parameters in the leptonic mixing matrix and in particular on 
the Majorana phases, a dependence that completely cancels out in vanilla leptogenesis. 
In this way flavour effects can potentially improve the predictive power of leptogenesis, 
though one still needs some parameter reduction as discussed in the unflavoured case. 
An extreme case occurs when the first term in Eq. vanishes for €, = 0, corresponding 

to total lepton number conservation in Nj,-decays. All the asymmetry has then to be 


generated by the second term [175]. In this case one can show easily that necessarily the 


38The tree level probabilities p?, ~ (Pia + Pig) /2: 
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Figure 9: Allowed region in the plane M, versus mı when flavour effects are included 
(from |181|). In the left panel no fine-tuning in the seesaw formula is allowed, in the right 


panel a mild fine-tuning is allowed. 


Apia s stem uniquely from the CP violating low energy neutrino phases, mainly from the 
Majorana phases, and that successful leptogenesis is attainable (158][182/183). 

For this reason, calculating the final asymmetry within a certain model including 
flavour effects is absolutely necessary. However, despite this, the absolute bounds on M: 
and my, we discussed within vanilla leptogenesis and shown in Fig. 8 are not dramatically 
relaxed in the hierarchical case. The reason is that in the absence of wash-out, as we 
mentioned, one recovers the unflavoured regime and, in the presence of wash out the 
relaxation is large but, unless some fine-tuning in the seesaw formula is introduced, is 
not at the level to relax the lower bound on M, considerably and similarly for the upper 
bound on my. However, if some fine-tuning in the seesaw formula is allowed then, even 
for hierarchical RH neutrinos, one can lower the scale of leptogenesis. In Fig. (9) we show 
a comparison between the lower bound on M, when flavour effects are included compared 
to vanilla leptogenesis. In the left panel no fine-tuning in the seesaw formula is allowed 
while in the right panel a mild fine-tuning is allowed and one can see how the lower bound 
is relaxed by one order of magnitude. This possibility of lowering the scale of leptogenesis 
has been also implemented in a model and it has been also recently thoroughly 
studied in (185). 

Notice that if Mı < 10° GeV interactions of leptons produced by RH neutrino decays 


with the RH muons in the thermal bath are effective in decohering the e + u component 
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of the |@;)’s as well. In this situation one has to consider a three fully flavoured regime 
where now one has to track the flavoured asymmetries in all three flavours. Notice that 
for values of Mı ~ 10'2 GeV and Mı ~ 10° GeV one is in a sort of transition between the 
unflavoured and the two fully flavoured regime and between the two fully flavoured regime 
and the three fully flavoured regime respectively. In these transition regimes the validity 
of Boltzmann equations breaks down and one needs to use a density matrix formalism to 
describe the transition (67\/176)/185}{187]. 

N2-leptogenesis. We have so far assumed that the asymmetry is dominated by 
the contribution of Nj-decays. As far as flavour effects are neglected, in the case of 
vanilla leptogenesis, this assumption works in most cases with just some special exception 
relying on having Kı < 1, in which case one can have that the asymmetry produced by 
next-to-lightest RH neutrinos can dominate (N2-leptogenesis) |188]. 

However, when flavour effects are taken into account, the asymmetry produced from 
the next-to-lightest (2) RH neutrinos can find its way to survive the N, wash-out 
much more easily and N2-leptogenesis cannot be regarded as a special case any more. 
Therefore, the contribution to the final asymmetry from Nə-decays needs in general to be 
systematically taken into account, unless special assumptions are made, like for example 
that Mə > Tro > Mı. 

There are two reasons why the asymmetry from N>-decays can survive more easily. 
The first reason is that for Mı < 10° GeV the lightest RH neutrino wash-out occurs in the 
fully three flavoured regime. At the lightest RH neutrino wash-out stage, the asymmetry 
produced from N2-decays can be regarded as a pre-existing initial asymmetry. However, 
now the Eq. holding in the unflavoured case, has to be written as the sum of three 
contributions, one from each flavour [180] 


f i 3r iSi on 
NE r = Nx. e 8 Kie + NK e 8 Kip + NR e 8 Kir . (147) 


In our case the three flavoured pre-existing asymmetries are produced from Nə-decays. 
Now, having included flavour, it is sufficient that just one Kia < 1 for the asymmetry 
in that flavour to survive. In the same random generation of models performed in [160], 
where it was found that the probability to have K, < 0.1 is approximately 0.1%, it was 
also found that the probability for Kie < 1 is 26% and the probabilities for Kip, Kır <1 
are ~ 7%. Therefore, it is clear that in this case the survival of an asymmetry produced 
from Nə-decays is not special but it becomes a very plausible situation. As we will discuss 
in the next subsection, this situation is crucial for SO(10)-inspired models. 

A second reason why the asymmetry from N -decays, or at least some fraction of it, can 


survive Nj-decays is simply due to flavour projection: considering the full three-dimensional 
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structure of flavour space, the wash-out from Nj, inverse processes acts only on a component 
that is parallel to the 4 flavour, i.e., the flavour of leptons produced in N,-decays and 
that can induce N; production from Nj, inverse decays. The orthogonal component will 
completely escape the wash-out (67|[189}. This effect is interesting since it holds even 
if the N, and Nə masses fall in the same flavour regime and, in particular, even in the 
unflavoured regime, for Mı, Ma >> 10!? GeV. An interesting application is that a new 
allowed region, where the No asymmetry dominates, opens up in the space of parameters 
in the two RH neutrino case {190}. 


5.4 SO(10)-inspired leptogenesis 


One of the most intriguing application of flavour effects is that they rescue SO(10)-inspired 
leptogenesis (191). Given the typical very hierarchical RH neutrino mass spectrum rising 
in SO(10)-inspired models, see Fig. |7| one has Mı < 10°GeV and, as we already 
pointed out, the asymmetry produced from the lightest RH neutrinos would be orders 
of magnitude below the observed value. Moreover, in SO(10)-inspired models one has 
Kı > 1. Considering the expression Eq. for the wash-out of an asymmetry 
generated prior to the lightest RH neutrino wash-out assumed to occur before sphalerons 
go out-of-equilibrium, the asymmetry produced from N2-decays would be efficiently washed 
out and N2-leptogenesis would not be viable. 

However, when flavour effects are taken into account, the expression Eq. for the 
final value of a pre-existing asymmetry gets modified into Eq. and, specifically in 


the case of N2-leptogenesis, when N2-decays take place in the two fully flavoured regime 


for 10°? GeV < Mz « 10! GeV, one has 


e Kae Koe pn 
NET A | 2 Ent K Ko, ) + (a = en TA) a 3 Kie d 
275 


K K 3T 
cis, | 2u Eart K(Kor1) + (a = = en ig 2) Bos Kiun Ji 


+ orka T (148) 


where Kor = Ko + Ko, and Eorl = E2e + Cay: In this expression we have also included 
so called phantom terms [193] that, however, in the case of SO(10)-inspired leptogenesis 
give a very small contribution and can be neglected. 

The interesting property is that, analogously to vanilla leptogenesis where the final 
asymmetry depends only on the high energy parameters associated to the lightest RH 
neutrino Nj, the final asymmetry in SO(10)-inspired leptogenesis depends only on the high 
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energy parameters associated to the next-to-lightest RH neutrino Nz. The flavoured CP 
asymmetries and the flavoured decay parameters can be expressed through the parameterisation 
we discussed in 4.4.2 (see Eqs. and (110)). Out of the three a;’s, the asymmetry 
depends only on a2. Moreover, since the matrix V; ~ Verm œ I, the dependence on the 

six parameters in the matrix Vz is mild. For this reason, from the successful leptogenesis 
condition, one is able to obtain predictions on the low energy neutrino parameters, and 

to place interesting experimental constraints, in particular [191}{193): 


e Inverted ordering is now in strong tension with the cosmological upper bound on 


neutrino masses (see (Eq. [84)); 
e There is a lower bound 13 Z 2°, now nicely confirmed by the 013 discovery; 
e There is a lower bound on the lightest neutrino mass mı 2 1 meV; 


e Majorana phases are constrained within narrow ranges around special values and 
for this reason the effective 0vGG mass Mee cannot be arbitrarily small and bulk 


solutions are found for values above meV. 


Recently it has been shown how the lower bound on the absolute neutrino mass 
scale, both on mı and Mee, becomes much more stringent for values of ô and 693 in a 
region (in orange in Fig. including current best fit values of ô and 623 (see Eq. (80)), 
and from Fig. [10] one can see that one obtains 


mı Z 34meV and Mee 2 3lmeV. (149) 


This is quite interesting since future experiments should be able to fully test these lower 
bounds. It is in particular interesting that, despite these solutions are found for normally 
ordered neutrino masses, there is such a stringent lower bound on Mee. These results 
clearly show that even high-scale leptogenesis models are testable. From this point of 


view the discovery of 0v36 decay would be a clear (necessary) breakthrough discovery 


supporting leptogenesis [°] 


39While on the other hand notice from Fig. |10| how a (low energy) CP conserving value 6 = 7 is 
perfectly compatible with successful leptogenesis. This is a clear example showing how observing CP 
violation in neutrino oscillation experiments is neither a necessary nor sufficient condition for successful 
leptogenesis. However, of course, one would in general expect to observe it, while ruling it out with 
great precision would represent quite a special situation casting some doubts, also because there are no 


compelling models where this is predicted to occur. 
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! m4 2 34 meV 


Figure 10: Lower bound on m; (left panel) and Mee (right panel) from SO(10)-inspired 
leptogenesis (from |194|). 


6 Dark matter from active-sterile neutrino mixing and 
the vMSM model 


As we discussed in the previous section, leptogenesis provides an elegant way to the 
solution of the matter-antimatter asymmetry of the universe. In its minimal, and original, 
version it just relies on the type-I seesaw extension of the SM Lagrangian Eq. that 
also successfully explains the observed neutrino masses and mixing. It is then legitimate to 
wonder whether the same extension might also address the dark matter of the universe. 
It is quite intriguing that there is indeed a very simple solution relying on the 
active-sterile neutrino mixing, expressed by Eq. (106), predicted by the type-I seesaw 
Lagrangian. 

The possibility for the heavy seesaw neutrinos to provide a viable dark matter candidate 
necessarily relies on special regions of parameter space where they are stable on cosmological 
scales. Their lifetime needs necessarily to be longer than the age of the Universe to be 
able to play the gravitational role of dark matter as cosmic glue, explaining the dynamics 
of clusters of galaxies and galactic rotation curves. However, depending on their mass, 
the lower bound on the life-time is orders of magnitude more stringent, in order to avoid 
experimental constraints deriving from indirect searches and also, at some energies, from 
cosmic microwave background and 21 cm global cosmological signal. 


If the mass of the heavy neutrino Ny is higher than the Higgs boson mass, then its 
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lifetime, taking the inverse of the decay width in Eq. (124), is simply given by 


TI = 


8m 5 (FF 


2 
~ 10 150 
(hth), M, Kı \ M; ) a rd ey) 


that we wrote in terms of the decay parameters defined in Eq. (16). Therefore, one can 


see that the stability condition Ty > tọ ~ 4 x 10! sec, would translate into an upper 


bound 


(151) 


In this limit, the RH neutrino would decouple from the seesaw formula and one would have 
simply mı = Mr x 1/M;. Within this minimal extension of the SM there cannot be any 
production mechanism relying on such tiny Yukawa couplings, even barring fine-tuning 
considerations on the feasibility of building such a model. One could resort to some 
gravitational production mechanism but clearly this would lack predictive power being 


disconnected from neutrino parameters. 


6.1 The Dodelson-Widrow mechanism 


There is, however, a second possibility to consider. One RH neutrino could be lighter than 
gauge bosons and in this case it would not directly decay through their Yukawa couplings 
into gauge bosons and leptons. However, it would still decay into SM particles because of 
the active-sterile neutrino mixing since, as expressed by Eq. (106), heavy seesaw neutrinos 
have a LH component. The dominant decay would be the decay into three leptons. The 
possibility to have 7r > to can be realised if Mr < 2m,. One would still need small Yukawa 
couplings with Ky < 1 but not as tiny as in Eq. (151). In this case, the dark matter 
heavy neutrino has to be necessarily Nj, the lightest of the three, in order for the two 
heavier ones to reproduce the solar and atmospheric neutrino mass scales from the seesaw 
formula [197]. The lightest seesaw neutrino mass would then be inversely proportional to 
the lightest ordinary neutrino mass mı. The general expression Eq. for the mixing 


can then be specialised to the case when one has a vı — N; dominant mixing: 


[Ufa + (Uavra)®| — [Via Sdn Ym + (Via Eta vm)" (152) 
Ni = [Yar t+) [fiarta + (iarra) ] - (153) 


Vy 


This mixing would be then responsible for its decay dominantly into three neutrinos with 


M. —5 1 —4\ 2 
Tı ~ 5 x 1078s (5) (>) ; (154) 
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a lifetime 


where 6 = X, |Mpo1/M,|? is an effective active-sterile neutrino mixing angle and where 
one can notice the strong inverse-fifth power dependence on the mass typical of three 
particle decays (Sargent’s law). The expression is written in a way to highlight that values 
6 ~ 1074 are necessary in order to have such long lifetimes that would escape stringent 
experimental constraints from X-ray observations in galaxy clusters. The striking coincidence 
is that this order-of-magnitude for the value of the mixing angle is also the correct one 
able to reproduce the observed dark matter abundance (see Eq. (47). Indeed, solving a 
simple Boltzmann equation for the production of the N,’s, one finds for its contribution 


to the energy density parameter 


0 i M i 
2 — BaS 
Qy,h 0.1 (Ga) (= ) F (155) 


that can be indeed satisfied for 0 ~ 1074 and Mı ~ keV. This is often referred to as the 
Dodelson-Widrow mechanism for the production of a sterile neutrino abundance able to 
reproduce the observed dark matter abundance [195]. Interestingly, for keV masses, the 
N,’s would behave as warm dark matter, implying that compared to cold dark matter 
there would be a reduced power in the large scale structure at scales corresponding to 
dwarf galaxies (~ 0.1 Mpc in comoving length). This would also help in solving some 
potential issues in pure cold dark matter N-body simulations that seem to predict too 
many satellite galaxies, in a galaxy like ours, compare to what is observed astronomically. 
At the same time warm dark matter would also smooth the cusp profile in galaxies that 
is predicted by cold dark matter N-body simulations but that is in tension with different 
observations. 

Sub-dominantly the N,’s can also radiatively decay and this offers a very beautiful 
experimental signature to test them with X-ray observations that place an upper bound 
M, < 10keV for 6 ~ 10~4 (see [198|/199}). 

Moreover large scale structure N-body simulations place a lower bound Mı = keV, 
since otherwise it would be too warm. These constraints have progressively ruled out 
the minimal scenario relying on the assumption of a non-resonant active-sterile neutrino 
mixing [199}[200], holding in the reasonable case that the initial lepton asymmetry is of 
the order of the baryon asymmetry (as for example it occurs in the minimal scenario of 
leptogenesis). However, if a large lepton asymmetry (L > 1075) is present at the time 
of mixing, then the mixing would be resonant and the production is enhanced for the 
same mixing angle compared to a non-resonant production (201]/202]. In this way one can 
reproduce the correct abundance and at the same time having longer life-times, reconciling 


predictions with the experimental constraints. 
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Interestingly, the observation of a new 3.5 keV line in the X-ray observations of different 
clusters of galaxies could be explained by such a mechanism for a decaying 7 keV 
sterile neutrino with a mixing angle 0 ~ 4 x 10~° and its relic abundance can explain 
the observed dark matter density parameter for a lepton asymmetry at the resonance 
L~4.6 x 10-4 [205]. Recently, observations from blank-sky observations testing emission 
from the dark matter Milky Way halo seems to have excluded the interpretation of the 
3.5 keV anomaly in terms of a 7 keV decaying sterile neutrino dark matter [206]. However, 
such a conclusion has very recently been questioned and more observations will be 


likely necessary to definitively test this exciting possibility. 


6.2 The vMSM model 


The non-resonant production of sterile neutrino dark matter can be nicely embedded 
within the seesaw mechanism in a way that neutrino masses and mixing parameters are 
also explained |197|. Even more ambitiously, one can also explain the matter-antimatter 
asymmetry of the universe via leptogenesis from the mixing of the two heavier RH 
neutrinos, the so-called ARS scenario (1721174), obtaining a unified model of neutrino 
masses and mixing, dark matter and baryogenesis: the so-called v MSM model [208]. 

The necessity of a resonant production with the generation of a large lepton asymmetry 
prior to the sterile neutrino production is clearly more challenging to reproduce but it 
can also be accommodated within the vMSM, assuming just a type-I seesaw extension of 
the SM with three RH neutrinos. This is possible if the two heavier RH neutrinos Nə and 
N; produce the large lepton asymmetry after sphaleron processes go out-of-equilibrium 
(by mixing or decays) and clearly prior to the lightest RH neutrino production. This 
late leptogenesis would then be instrumental in the sterile neutrino dark matter resonant 
production [209}211]. 

The possibility to generate a large lepton asymmetry prior to the dark matter imposes 
very stringent constraints. First of all the two heavier RH neutrinos have to be extremely 
degenerate with a mass difference that has to be in an extremely narrow range of values 
around M — Mz, ~ 107! Mp, a condition rising from a necessary strong fine-tuning 
between light and heavy neutrino mass differences [212]. Moreover, in order to have 
sufficient CP violation one also obtains a lower bound (holding for normal ordering) 
Mə 2 2 GeV. Despite the strong necessary fine-tuning, this scenario is quite interesting 
since it can be tested either from the decays of the dark matter particles that, again, 
should produce a line in X-ray observations, and/or also from B meson decays, that might 


produce the two heavier RH neutrinos if these are lighter than ~ 5 GeV, in experiments 
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such as SHiP [213], MATHUSLA and the approved FASER [215]. Clearly the 
observation of the 3.5 keV line has triggered great excitement. It should be said however 
that a recent numerical analysis where the mass of Nı has been fixed to 7keV, the 
value that would explain the anomalous 3.5 keV line, has found no solutions in the 
space of parameters able to reproduce both the baryon asymmetry and the observed dark 
matter energy density parameter [216]. The main problem is the difficulty to produce, 
after sphaleron out-of-equilibrium and prior to dark matter production, the large lepton 
asymmetry needed to enhance resonantly the dark matter production in active-sterile 
neutrino mixing. In order to achieve dark matter resonant production, solutions where 
such large asymmetry is produced from the N23 late decays have been recently found 
but an exquisite degree of fine-tuning of the involved parameters is required. Indeed, 
the necessary lepton asymmetry for the resonant production is generated if the physical 
splitting (M3 — M2) ~ 10716 M, with Mz ~ GeV and if the real part of the complex angle 
in the orthogonal matrix is close to +7/2 with a relative precision of ~ 1078 for normal 
hierarchy and ~ 107” for inverted hierarchy (see Fig. 2 in (217]). These solutions are also 
compatible with successful leptogenesis from RH neutrino mixing. If such fine-tunings are 
not imposed, then a generation of just (5-10)% of the observed dark matter abundance 


is generically obtained for Mə 3 ~ GeV. 


7 Beyond the minimal seesaw Lagrangian 


The difficulty, if not impossibility, to obtain a scenario on the origin of matter in the 
universe within a minimal type-I seesaw extension of the SM, motivates an extension 
of the scenario. Moreover, unless an origin of the 3.5 keV anomaly in terms of sterile 
neutrino decays is confirmed | the strong existing constraints from X-ray observations 
and structure formation seem to disfavour a solution of the dark matter puzzle in terms 
of a ~ keV seesaw neutrino. This encourages the investigation of an extension of the SM 


beyond the minimal type-I seesaw. We briefly discuss a few recent proposals. 


7.1 Higgs induced RH neutrino mixing (RHiNo) model 


The success of minimal leptogenesis in explaining the matter-antimatter asymmetry of 
the universe in relation to the observed neutrino masses and mixing, motivates to retain 
the type-I seesaw Lagrangian as a successful starting point and considering an extension 
to it. 

40The launch in 2022 of the XRSM satellite should definitively solve the controversial situation (218). 
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Even more, one can think to repeat, in an upgraded version, the same successful 
strategy. We saw that the type-I seesaw is the simplest way to realise a UV-completion of 
the SM that corresponds to a low energy effective theory described by the 5-dim Weinberg 
operator (115). Having now introduced RH neutrinos, leading to an explanation of 
neutrino masses and mixing via seesaw mechanism and matter-antimatter asymmetry via 
leptogenesis, one can form a similar new 5-dim operator, the Anisimov operator 
(I,J = 1,2,3) 

04 =“ of ON Ny, (156) 
a variant of the Higgs portal operator. Having introduced this operator, one can go 
back to the idea of having one heavy RH neutrino, with vanishing Yukawa coupling as 
a candidate of cold dark matter. The Anisimov operator interestingly would induce, via 
the Higgs, a mixing between a seesaw heavy neutrino and the dark matter RH neutrino. 
This mixing would be responsible both for the dark matter RH neutrino production 
at high temperatures but also for its decays at zero temperature. This is the Higgs 
induced RH-RH neutrino mixing (RHiNo) model. The interesting feature of the model 
is that at high temperatures the mixing is enhanced œ T° by thermal effects. In this 
way the effect of the Anisimov operator is such to have an enhanced production at high 
temperatures but would fade out today, still however inducing dark matter decays, with 
the possibility to test them in indirect searches. Clearly the possibility to satisfy stability 
while still reproducing the observed dark matter energy density parameter has to be 
verified quantitatively and goes through an accurate calculation both of the dark matter 
production and of all decay channels contributing to the total decay width of the dark 
matter RH neutrino and to indirect search signals |220]. 

Initially, calculations of the final relic dark matter abundance have been done imposing 
a resonance condition, leading to the condition that the dark matter is heavier than the 
source RH neutrino, and then using a Landau-Zener formula for the description of the 
non-adiabatic production of dark matter RH neutrinos at the resonance. However, recent 
results within a density matrix formalism show that the Landau-Zener formula 
fails to reproduce the correct dark matter abundance unless there is a strong degeneracy 
between dark and source RH neutrino. The problem is that the Landau-Zener assumes 
that oscillations have already developed at the resonance. However, this is true only if 
AM/Mpm « 107? Mpm/GeV, where AM is the mass difference between the dark matter 
and the source RH neutrinos. This condition clearly requires the two RH neutrinos to be 
quasi-degenerate. For hierarchical RH neutrinos, density matrix calculations show that a 


non-resonant production is still occurring but the final abundance is order of magnitudes 
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lower than within a Landau-Zener approach. This non-resonant production has analogies 
with the extremely non-adiabatic limit in solar neutrinos. In this case essentially thermal 
effects generated by the Anisimov operator are still important but those generated by the 
Yukawa coupling of the source RH neutrino, corresponding to medium effects in the case 
of solar neutrinos, are negligible. 

The results is that for Mpm > Ms > Muiggs ~ 125 GeV one can reproduce the correct 
measured dark matter energy density parameter only assuming the thermalisation of 
the source RH neutrino at the time of the dark matter production from the mixing. 
This thermalisation however cannot be achieved by Yukawa couplings so somehow one 
needs to resort to some external mechanism. However, for dark matter masses below 
the Higgs and gauge boson masses a two body decay channel, Npm > Ns > A+ 4, 
is kinematically forbidden and this allows to have larger source RH neutrino Yukawa 
couplings still satisfying dark matter stability on cosmological scales. In this case however, 
the usual scenario of leptogenesis from decays cannot be successful for the generation of the 
matter-antimatter asymmetry and one needs to resort to the ARS scenario of leptogenesis 
from mixing [172], as in the yMSM model. 

On the other hand, since the production is non-resonant anyway, a condition Mpm > 
Mg is not needed any more and one can explore the case Mg > Mpm. This might open a 
new interesting solution where the dark matter cam be the lightest RH neutrino species 
N, and the two heavier ones, Nj and N3, can produce the matter-antimatter asymmetry of 
the Universe just within the standard leptogenesis scenario and generate neutrino masses 
in two RH neutrino seesaw models. 

It is interesting that the mechanism can also be combined with a phase transition 
driven by a scalar that is also responsible for the generation of the source RH neutrino 


mass. During the phase transition both dark matter, neutrino masses and gravitational 


waves are generated [222]. 


7.2 Frogatt-Nielsen model and dynamical Yukawa couplings 


The idea that a heavy RH neutrino can play the role of dark matter and be produced 
through couplings that were much stronger at high temperatures but almost vanish at 
present time to ensure stability, can be exported to other models. Recently, it was 
proposed that Yukawa couplings of dark matter RH neutrino in the type-I seesaw Lagrangian 
could dramatically change during a phase transition of a new scalar X in a Frogatt-Nielsen 
model, passing from O(1) couplings to tiny values consistent with cosmological stability 


223|. Before the / phase transition the large Yukawa couplings thermalise the dark 
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matter RH neutrino. At the phase transition the dark matter RH neutrino Yukawa 
coupling is strongly suppressed and this induces the freeze-out of the dark matter RH 
neutrino abundance that has to be such that to yield the observed dark matter energy 
density parameters. In addition to the dark matter RH neutrino, the presence of two 
additional RH neutrinos can reproduce the neutrino masses and mixing via type-I seesaw 


mechanism and the matter-antimatter asymmetry via leptogenesis. 


7.3 Seesaw and dark side sector 


The RHiNo model is the first example of a model connecting neutrino masses 
from type-I seesaw mechanism to dark matter. Since it introduces a non-renormalisable 
operator, it is not UV-completed. Notice that the dark sector is represented by the RH 
neutrino Np». 

Along similar lines one can connect the seesaw physics to dark matter introducing a 
dark sector made of a new singlet complex scalar ¢@ and a singlet fermion x and imposing 
a Z4 symmetry [224]. The SM Lagrangian is then extended by 


1_~ fe 
L—-Lu = 7 LH ON + 5 MNON +ydXx+V(H, 4). (157) 


Prior to the Z4 symmetry breaking the dark matter remains massless and does not couple 
to the SM but after symmetry breaking it can be shown that ø mixes with the SM Higgs 
and dark matter can be produced and annihilate through processes SM + SM & x +X. 
The abundance can then be calculated via standard freeze-out solution. The dark matter 
is stable due to a relic Z symmetry. 

With the same dark sector a variation of the RHiNo models that allows to obtain a 


UV-completed model can be obtained introducing a neutrino portal term |225 
Lportal = YDS OX Nr T h.c. (158) 


In this case there is no mixing of @ with the standard model Higgs and the dark matter 
production occurs via freeze-in solution. In both models the DM particle is stable so 
that there are no indirect search signatures and this can be seen as a phenomenological 


drawback. 


7.4 Scotogenic model 


In the scotogenic model the SM is extended by adding three copies of SM singlet 


fermions plus an additional Higgs doublet that is inert since it does not develop a vacuum 
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expectation value. The new sector is odd under an unbroken Z symmetry so that the 
lightest Zə odd particle, that turns out the lightest neutral scalar of the inert doublet, can 
play the role of a stable dark matter particle candidate. The relic abundance is obtained 
via the standard freeze-out mechanism. The Z symmetry and the corresponding charges 
of the fields prevent the inert Higgs doublet to couple with SM fermions at tree level so 
that light neutrino masses are forbidden at tree level (otherwise one would have a seesaw 
generation) but can be generated radiatively at one loop. 

Interestingly, the out-of-equilibrium decays of the three SM singlet Zə odd fermions 
can also generate the baryon asymmetry via leptogenesis [227]. As in the case of vanilla 
leptogenesis, the asymmetry would be dominated by the decays of the lightest Zə odd 
fermion N, with mass M, and generated at temperatures T ~ Mı. It is found that 


successful leptogenesis can be obtained for Mı ~ 10 TeV. 


8 Final remarks 


Understanding the origin of matter in the universe poses a challenge to our current 
knowledge of fundamental physics, since within our current description of fundamental 
interactions based on general relativity and the SM there is no model able to account 
for the dark matter and matter-antimatter asymmetry of the universe. This is extremely 
interesting since a model of the origin of matter in the universe necessarily implies the 
existence of new physics. We mainly discussed solutions relying on an extension of the 
SM and we have particularly focussed on those solutions that can also address neutrino 
masses and mixing, since in the absence of clear new physics from other phenomenologies, 
this is an attractive path. From this point of view, the continuation of the already well 
established programme of low energy neutrino experiments might provide an important 
insight in future, especially in connection with leptogenesis. In this respect, the discovery 
of Ov decays would represent a clear breakthrough, establishing the Majorana nature 
of neutrinos and lepton violation at tree level. At the same time, the measurement of the 
Ov 66 effective neutrino mass would provide a way to test specific models and we have seen 
in particular the case of SO(10)-inspired models. When 0v8 experiments are combined 
together with the cosmological information on the absolute neutrino mass, one obtains a 
powerful test that might even lead, in optimistic scenarios, to pin down specific scenarios 
of leptogenesis. We have also seen how the 3.5 keV anomaly would in case point also to 
a solution of the dark matter puzzle based on neutrino physics. On this front, planned 
laboratory experiments, such as SHiP, MATHUSLA and FASER will test the existence of 


RH neutrinos lighter than ~ 5 GeV that are predicted in scenarios on the origin of matter 
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such as the vMSM. However, in the last years a wide variety of different experimental 
discoveries open new ways to investigate models of the origin of matter in the universe. 
The discovery of very high energy neutrinos at IceCube provides a new investigative tool 
on scales much above the TeV scale, unaccessible to ground laboratory experiments. The 
first experimental information on 21 cm cosmological global signal allows already to place 
new constraints on dark matter models probing a so far unexplored stage in the Universe 
history, the so-called dark ages at redshifts z = 10. Finally, the discovery of gravitational 
waves opens an incredible variety of new tests on models associated to the production of 
stochastic backgrounds and/or primordial black holes. 

Clearly history teaches us that one should not be overly optimistic about getting close 
to the solution to the problem of the origin of matter in the universe: too many times great 
announcements have then been followed by great disappointments. However, it is fair to 
say that we are certainly entering a new era where the variety of theoretical models, often 
seen as a sign of a perilous navigation in uncharted waters, is at least counterbalanced 
by an incredible variety of new phenomenological tools that can guide us through the 
journey, making everything extremely exciting. For these reasons, finding a solution to 
the origin of matter in the universe should certainly inspire new explorers eager to embark 


in such a journey toward new lands. 
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